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PHILIPS BUSINESS COMMUNICATIONS 


Internet Protocol (IP) Telefonie is een technologie waarmee spraakcommunicatie over 
datanetwerken, zoals het Internet, mogelijk wordt gemaakt. De belangstelling voor IP- 
Telefonie is enorm. Dit komt in de eerste plaats door de grote besparingen op telecom- 
municatiekosten die men verwacht te kunnen behalen. Op iets langere termijn zullen 
het vooral de toegevoegde diensten zijn die de grootste drijfveer gaan vormen voor IP- 
Telefonie. In dit artikel wordt stilgestaan bij de vraag hoe IP-Telefonie kan worden 
geïntegreerd in de bedrijfsomgeving waar van oudsher gebruik wordt gemaakt van 
PBX en PBX-netwerken. Eerst echter zal enige aandacht worden gegeven aan het fe- 
nomeen van IP-Telefonie. Wat is het, hoe werkt het en wat zijn de standaarden. 


Convergentie van IP-Telefonie en Inleiding 
systemen voor bedrijfs- 


telecommunicatie Het Internet brengt enorme veranderingen teweeg in 
de manier waarop we werken en communiceren. 
5 Mensen over de hele wereld hebben toegang tot het 
Internet en kunnen onderling informatie uitwisselen. 
Hierbij wordt gebruik gemaakt van diensten zoals 
elektronische post, file transfer en het World Wide 
Web. Video en audio kunnen probleemloos over het 
Internet worden getransporteerd en het aanbod van 
multimedia diensten op het Internet neemt snel toe. 


Internet Telefonie, of beter, Internet Protocol (IP) 
Telefonie is een technologie waarmee spraakcom- 
municatie vanaf de multimedia PC mogelijk wordt. 
Het geeft de gebruiker de mogelijkheid om over de- 
zelfde Internetverbinding tegelijkertijd telefoon- 
gesprekken te voeren en andere Internet-applicaties 


IN DIT ARTIKEL te gebruiken, bijvoorbeeld een Web-browser. 
e Internet, Intranetten en het Internet Protocol IP-Telefonie is niet beperkt tot PC-gebruikers. Met 
e __IP-Telefonie en IP-Telefonie Gateways behulp van IP-Telefonie Gateways is het mogelijk 


om de traditionele circuitgeschakelde telecom- 
municatienetwerken te koppelen met _pakket- 
gebaseerde IP-netwerken. Daarmee werpt zich ook 


e Netwerken voor bedrijfscommunicatie 
e _ Toepassingen van IP-Telefonie in de bedrijfs- 


omgevin, d 4 5 

Ë EE Virtual Networking over IP de vraag op, hoe IP-Telefonietechnologie geïnte- 
a Telewerken greerd kan worden in systemen voor bedrijfscommu- 
< Web Call Centers nicatie voor het bieden van nog betere diensten aan 


e De evolutie van de PBX de zakelijke gebruiker. 


1 Tijdschrift van het Nederlands Elektronica- en Radiogenootschap deel 64-nr.1-1999 


IP-Telefonie achtergronden 


Internet, intranetten en het Internet Protocol 

Het Internet is een wereldwijde verzameling van 
computernetwerken die met elkaar verbonden zijn 
door gebruik te maken van de Internet Protocol suite. 
Onderdeel van de IP-standaarden is een unieke ma- 
nier om alle deelnemende computers te adresseren 
met een naam. De belangrijkste karakteristiek van 
het Internet is dat het zo wijdverspreid is. Mensen 
over de hele wereld hebben toegang tot het Internet. 
De toegang wordt aangeboden door de Internet Ser- 
vice Providers. 


Intranetten werken op dezelfde manier als het Inter- 
net. Waar het Internet toegankelijk is voor het grote 
publiek, zijn intranetten voor exclusief gebruik door 
besloten gebruikersgroepen. Het voordeel van intra- 
netten is dat de capaciteit en de prestaties beter con- 
troleerbaar zijn. Dat komt doordat intranetten ge- 
woonlijk in eigen beheer zijn van de gebruikende 
organisatie. 


Zowel het Internet als de intranetten zijn pakketgeba- 
seerde netwerken, waarbij het Internet Protocol ver- 
antwoordelijk is voor het routeren van de informatie- 
pakketten door het netwerk. 


Protocol Distributie in Wide Area Networks 
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TL IP Distributie 
60% 
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‘Bron Garmer Groep 
Figuur 1: Het belang van het Internet Protocol 


Het belang van IP wordt duidelijk als we kijken naar 
figuur 1. Hieruit blijkt dat in het jaar 2000 wereld- 
wijd 80% van alle datanetwerken gebaseerd zal zijn 
op het gebruik van IP. 


IP-Telefonie 

IP-Telefonie is een technologie die het mogelijk 
maakt om telefoongesprekken te voeren over het In- 
ternet, een intranet of elk ander op IP-gebaseerd 
netwerk. In tegenstelling tot de traditionele circuitge- 
schakelde telefonie, maakt IP-Telefonie gebruik van 


De Internet Protocol (EP) suite 6 | 


Applicatieprotocollen Applicatieprotocollen 
bijv. HTTP, FTP, TELNET bijv. SNMP, TFTP | 


Transmission Control Pro- 
tocol (TCP) voor betrouw- 
baar transport 


Ë 
| 
Í 


User Datagram Protocol 
(UDP) voor 
onbetrouwbaar transport 


Internet Protocol (IP) 


Fysieke Netwerk 
(LAN, WAN) 


Terzijde | 


pakketgebaseerde technieken. Hierbij wordt de 
spraakinformatie opgedeeld in IP-pakketten welke 
over een IP-netwerk worden verzonden zonder voor- 
af een verbindingsweg op te zetten. De ontvanger 
verzamelt de IP-pakketten en stelt daarmee het oor- 
spronkelijke spraaksignaal weer samen. Het gevolg 
van deze methode is dat IP-Telefonie geen garanties 
biedt op voldoende en constante bandbreedte voor 
een gesprek. Dit kan resulteren in vertragingen en 
mindere kwaliteit van de spraakoverdracht. Tot op 
heden hebben alle IP-Telefonieproducten hinder van 
dit probleem. De verwachting is echter dat de kwali- 
teit snel en drastisch zal toenemen met nieuw be- 
schikbaarkomende producten. Verbeterde compres- 
sietechnieken zullen de hoeveelheid vereiste band- 
breedte voor een gesprek nog verder omlaag bren- 
gen. Aan de andere kant zal ook de hoeveelheid be- 
schikbare bandbreedte gaat toenemen door verder- 
gaande investeringen in de IP-netwerken. 


IP-Telefonie vindt zijn oorsprong in de consumen- 
tenmarkt. De multimedia PC van de gebruiker is uit- 
gerust met een microfoon, luidsprekers, een geluids- 
kaart, een speciaal softwarepakket, en natuurlijk een 
Internet aansluiting. 


Internet Service 
Provider 


Figuur 2: IP-Telefonie in de consumentenmarkt 


Het speciale softwareprogramma in de PC van de 
oproeper comprimeert en converteert de spraaksig- 
nalen in datapakketten. Deze datapakketten worden 
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vervolgens als IP-berichten via het IP-netwerk naar 
de PC van de bestemming gerouteerd. Bij de be- 
stemming worden de datapakketten weer terug ge- 
converteerd en geëxpandeerd tot de oorspronkelijke 
spraaksignalen. Het routeren van de IP-berichten 
door het netwerk is gebaseerd op het gebruik van 
IP-adressen (elke PC die gekoppeld is met het Inter- 
net krijgt een uniek IP-adres of een naam die ver- 
taalbaar is in een uniek IP-adres). 


De belangrijkste voordelen van IP-Telefonie zijn: 

e _ Kostenbesparing, onder andere door het omzei- 
len van gesprekskosten. 

e Toegevoegde waarde door simultaan gebruik 
van spraak, video en data en door een rijkere 
gebruikerservaring. 

e _ Gemak door gemeenschappelijke infrastructuur 
en randapparatuur voor spraak en data. 


Momenteel zijn de meeste toepassingen van IP- 
Telefonie nog beperkt tot het leveren van basistele- 
foniediensten, maar toegevoegde functies en 1we 
toepassingen worden snel belangrijk. Daarbij is IP- 
Telefonie zeker niet beperkt tot alleen maar het 
transporteren van spraak. 


Spraak over IP. Facsimile over IP 


Private Networking Audio conferentie 


Web Call Centres Telewerken 


Unified Messaging _ Video conferentie 


Figuur 3: Opkomende IP-Telefonie Markten 


IP-Telefonie technologie is ook geschikt voor het 
transport van facsimile over IP netwerken, voor ge- 
bruik in telewerkomgevingen, voor toepassingen in 
bedrijfscommunicatie netwerken en voor vele soor- 
ten van multimedia toepassingen. 


Standaarden voor IP-Telefonie 

De eerste producten voor IP-Telefonie waren ge- 
sloten producten. Dat wil zeggen dat twee gebrui- 
kers alleen met elkaar kunnen communiceren wan- 
neer ze op elkaar afgestemde producten van dezelf- 
de fabrikant gebruiken. Het was meteen al duidelijk 
dat dergelijke gesloten producten nooit kunnen lei- 
den tot een wereldwijde acceptatie van IP-Telefonie. 
Standaarden zijn daarom hard nodig om IP- 
Telefonie tot een succes te maken. Internationale 
standaardisatie is inmiddels dan ook in volle gang 
en eerste standaarden zijn reeds beschikbaar. De 
belangrijkste basis voor deze standaarden zijn de 
H.323 aanbevelingen van het ITU-T (zie de hier- 
naaststaande achtergrondinformatie). 


ITU-T aanbeveling H.323 


Systeem | Verbinding | Gatekeeper Audio/ 
besturings | besturings | (RAS) RTP spraak 
functies functies functies 6T.xx 
H245 H2250 H2250 | H225o | RTP 

H.225,0 
TP UDP 
Internet Protocol (IP) 


H.323 is een aanbeveling van het ITU-T. H.323 wordt gezien 
als de belangrijkste basis voor standaarden ten behoeve van 
multimedia communicatie (audio, video, data) over pakketgeba- 
seerde netwerken (bijv. een IP-netwerk) die niet noodzakelij 
kerwijs kwaliteitsgaranties bieden. H.323 is in feite een ver- 
zameling van een aantal andere aanbevelingen, zoals weergege- 
ven in de bovenstaande figuur. 


De belangrijkste componenten in het H.323 bouwwerk zijn: 


«__ _Systeembesturingsfuncties, vastgelegd in aanbeveling 
H.245. Bijvoorbeeld de onderhandeling tussen twee com- 
municerende eindpunten over het gebruik van een ge- 
meenschappelijke codec. 

«__Verbindingsbesturingsfuncties, vastgelegd in aanbeveling 
H.225.0. Deze omvatten de signaleringsfuncties voor het 
opzetten, onderhouden en verbreking van “gesprekken”. 


«__Gatekeeper functies, vastgelegd in aanbeveling H.225.0. 
De Gatekeeper is een optioneel onderdeel van een IP- 
Telefonienetwerk. Taken van een Gatekeeper zijn bij- 
voorbeeld adresvertaling en registratie en toelating van 
H.323 terminals in het netwerk (Registration, Admission, 
Status). 

« Real-time Transmission Protocol (RTP) en Realtime 
Transmission Control Protocol (RTCP), eveneens vastge- 
legd in aanbeveling H.225.0. Deze bieden functies voor 
het transport van real-time informatie over IP- 
netwerken. 


«___Audio/spraak functies. Voordat de audio informatie tus- 
sen de gebruikers kan worden getransporteerd, zal de in- 
formatie moeten worden gedigitaliseerd en gecodeerd. 
Hierbij kan gebruik worden gemaakt van verschillende co- 
decs (beide eindpunten in de communicatie dienen vooraf 
te onderhandelen over welke codec gebruikt gaat worden). 


«Video functies (optioneel) voor het transporteren van 
video informatie tussen eindgebruikers. 


H.323 componenten maken gebruik van de transportfuncties 
van hetzij TCP hetzij UDP, afhankelijk van de eisen die worden 
gesteld aan prestatie en betrouwbaarheid van het transport. 
Besturingsinformatie vraagt een betrouwbaar transport en 
maakt gebruik van TCP. Spraak/audio maakt gebruik van UDP 
omdat TCP veel complexer is en een nadelig effect heeft op de 
prestatie van de real-time communicatie. 


Terzijde 
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IP-Telefonie Gateway 

Om gesprekken te kunnen voeren tussen een Internet 
PC en de traditionele telefonie-omgeving is een IP- 
Telefonie Gateway nodig, zie figuur 4. 


Figuur 4: IP-Telefonie in de consumentenmarkt 


De belangrijkste taken van een IP-Telefonie Gateway 
zijn: 

« _ Het bieden van een interface naar het IP-netwerk. 
e Het bieden van een interface naar het telefonie- 
netwerk, bijvoorbeeld een PBX of het ISDN. 
Conversie tussen telefoonnummers en _ IP- 
adressen. 

Verzenden, ontvangen en converteren van spraak 
in IP-datapakketten. 


Bedrijfscommunicatienetten 


Bedrijven, organisaties en instellingen bezitten (of 
huren) vaak een eigen netwerk voor de interne tele- 
communicatie. Deze gesloten netwerken zijn daarbij 
gewoonlijk gekoppeld met openbare netwerken, voor 
de communicatie met externe gebruikers. Van ouds- 
her onderkent men twee soorten van bedrijfsnetwer- 
ken, een bedrijfstelefonienetwerk en een bedrijfsda- 
tanetwerk. Redenen voor het voeren van eigen be- 
drijfsnetwerken zijn dezelfde voor telefonie en data- 
netwerken: 
e Beheersing van de kosten (bij een eigen netwerk 
is de interne communicatie “gratis”). 
e Sneller en goedkoper systeembeheer door onaf- 
hankelijkheid van openbare dienstverleners. 
« Meer en op maat gesneden toegevoegde functies. 
e _ Betrouwbaarheid en veiligheid. 


De belangrijkste component in een bedrijfsnetwerk 
voor telefonie is van oudsher de PBX (Private 
Branch eXchange). De PBX maakt de verbindingen 
tussen de zakelijke gebruikers voor hun dagelijkse 
communicatie. Daarbij biedt de PBX een scala van 
toegevoegde functies die veel verder gaan dan wat 
het openbare telefonienetwerk kan bieden. Deze 
diensten zijn niet alleen beschikbaar voor gebruikers 


van vaste toestellen, maar ook voor mobiele gebrui- 
kers. Bijvoorbeeld door het gebruik van DECT 
(Digital Enhanced Cordless Telecommunications) als 
toegangstechniek tot de PBX. Voor gesprekken met 
de buitenwereld heeft de PBX koppelingen met het 
openbare telefonienetwerk. Daarnaast is het mogelijk 
om meerdere PBX’en te koppelen tot een groter be- 
drijfsnetwerk, bijvoorbeeld voor organisaties die uit 
verschillende vestigingen bestaan. 


Openbare netwerk 
(PSTN, ISDN) 


Figuur 5: Gekoppelde PBX'en in een bedrijfsnetwerk 


PBX'en in een netwerk kunnen zijn gekoppeld met 
huurlijnen of met geschakelde lijnen over het open- 
bare netwerk. Eén van de voordelen van het gebruik 
van huurlijnen is dat rechtstreekse signalering moge- 
lijk is tussen de PBX'en zonder interventie van het 
openbare netwerk. Dit opent de weg voor een hoge 
functionaliteit binnen het totale bedrijfsnetwerk. Het 
gebruik van huurlijnen is niet in alle configuraties 
economisch te verantwoorden. Het alternatief is dan 
om de PBX’en te koppelen via publieke geschakelde 
lijnen. Echter daarbij is er in eerste instantie geen 
rechtstreekse signalering mogelijk tussen de PBX'en. 
Daardoor is er ook geen transparantie voor PBX- 
diensten. Om dit probleem op te lossen, werd ge- 
zocht naar middelen om toch rechtstreekse signale- 
ring te kunnen voeren tussen PBX'en die zijn gekop- 
pelde door publieke geschakelde lijnen. De Philips 
SOPHO PBX bijvoorbeeld, biedt een “Private Virtu- 
al Networking” oplossing, waarbij PBX-signalering 
transparant kan worden getransporteerd over het pu- 
blieke netwerk. 


Figuur 6: Gekoppelde LAN's in een bedrijfs- 
datanetwerk 


Het hart van een bedrijfsnetwerk voor data is het 
Local Area Netwerk (LAN). LAN’s worden vaak 
toegepast om gemeenschappelijke hulpmiddelen 
(“servers”) zoals printers, mainframes en disks te 
delen tussen meerdere gebruikers. LAN's kunnen 
onderling worden gekoppeld door huurlijnen of indi- 
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rect via Wide Area Netwerken. Veel van de huidige 
LAN's maken gebruik van Ethernet verbindingen 
voor de fysieke koppeling naar de gebruiker. ATM, 
ook vaak genoemd als oplossing naar de desktop, 
lijkt daardoor vooral beperkt te gaan worden tot toe- 
passing in het Wide Area netwerk. Een belangrijke 
verdere constatering is, dat IP zich inmiddels stevig 
heeft gevestigd als het meest belangrijke netwerk- 
protocol voor communicatie in datanetwerken. 


Internet Protocol (IP) 


[ Ethernet ATM Frame Relay 


Figuur 7: IP boven op Ethernet, ATM, enz. 


Traditioneel zijn de bedrijfsnetwerken voor telefonie 
en data gescheiden netwerken. Niet alleen fysiek, 
maar ook organisatorisch. Waar de telecommuni 
tiemanager verantwoordelijk is voor het telefonie- 
netwerk, is de IT-manager verantwoordelijk voor het 
datanetwerk. Meer recent zien we echter dat het ver- 
schil tussen telefonie- en datanetwerken aan het ver- 
dwijnen is. Een voorbeeld hiervan is de opkomst van 
Computer Telephony Integration (CTI), een ontwik- 
keling waarbij een computer in het datanetwerk de 
gespreksafhandeling in de PBX kan beïnvloeden en 
besturen. En met een nieuwe techniek als IP- 
Telefonie wordt een belangrijke volgende stap gezet 
naar een verdere integratie van spraak- en datanet- 
werken. 


Toepassingen van IP-Telefonie in de 
bedrijfsomgeving 


De belangrijkste argumenten voor toepassingen van 

IP-Telefonie in de bedrijfsomgeving zijn de volgen- 

de: 

© Kostenbesparing op gesprekken en infrastructuur. 

« Toegevoegde waarde door gelijktijdige spraak, 
video en data en door een rijkere gebruikersper- 
ceptie. 

e Gemak voor beheerders en gebruikers door één 
enkele infrastructuur en één enkel gebruikersap- 
paraat (de PC). 


Bij het definiëren van IP-Telefonietoepassingen in de 
bedrijfsomgeving is het belangrijk om een volledige 
integratie na te streven met de bestaande telecommu- 
nicatie- en datacommunicatie-infrastructuren en sys- 
temen. Daarbij is het belangrijk dat gebruikers van 
IP-Telefonie kunnen beschikken over dezelfde dien- 


sten en faciliteiten als worden geboden in bestaande 
bedrijfstelefonienetwerken. 

In het nu volgende deel zullen enkele voorbeelden 
worden besproken van toepassingen van IP-Telefonie 
in de bedrijfsomgeving. De voorbeelden zijn afkom- 
stig van enkele succesvolle proeven die zijn opgezet 
met de Philips SOPHO PBX. Deze proeven passen in 
de plannen van Philips Business Communications 
om IP te voeren als het kernprotocol voor zijn be- 
drijfscommunicatieplatform. 


Private Virtual Networking over IP netwer- 
ken 

In deze toepassing worden de verschillende PBX'en 
van een gedistribueerd bedrijfstelefonienetwerk met 
elkaar gekoppeld via een IP-netwerk. Het is een uit- 
breiding van de al bestaande Private Virtual Net- 
working oplossing voor circuitgeschakelde verbin- 
dingen. De voordelen van deze uitbreiding zijn toe- 
gespitst op kostenbesparing op de dagelijkse be- 
drijfscommunicatie. Dit wordt bereikt door het ver- 
keer tussen de PBX'en te routeren over de bestaande 
datalijnen van het bedrijf (opvullen van de overcapa- 
citeit) of door het verkeer te routeren over het (nage- 
noeg kostenloze) publieke Internet. Hierdoor zijn 
minder huurlijnen nodig en minder kosten die moe- 
ten worden betaald aan de publieke operators. 


SOPHO PBX 
IP-Telefonie IP-Teleforie 
Goteway 


verbindingen} 
(ijv- ISDN } 
of huurlijn, } 
of Private 
Virtual} 
Networking à 
over ISDN) 


Figuur 8: Private Virtual Networking over IP 


De koppeling van de PBX met het IP-netwerk is door 
middel van een IP-Telefonie Gateway. Er zijn twee 
manieren waarop dit kan worden uitgevoerd. In de 
basisoplossing wordt alleen de spraak tussen de 
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PBX’en over het IP-netwerk vervoerd. Er is dan geen 
intelligente communicatie mogelijk tussen de 
PBX'en en daarmee is er geen transparantie van 
PBX-diensten. Beter is het om behalve de spraak ook 
de inter-PBX signalering te transporteren waardoor 
volledige transparantie van faciliteiten wordt verkre- 
gen, op dezelfde manier als wanneer gebruik zou 
worden gemaakt van huurlijnen. Dit is geïllustreerd 
in figuur 8. 

In deze toepassing is de gebruiker zich er niet van 
bewust als zijn gesprek over het IP-netwerk wordt 
gerouteerd. De handelingen die hij of zij moet ver- 
richten voor het opzetten van een gesprek zijn als 
gebruikelijk. De SOPHO PBX zorgt ervoor dat het 
gesprek, ongemerkt door de gebruiker, over de meest 
efficiënte verbinding wordt gerouteerd. Dat kan een 
IP-verbinding zijn, maar ook een huurlijn of een ver- 
binding over het geschakelde openbare telefoonnet- 
werk. De selectie van de meest geschikte route wordt 
voor elk gesprek opnieuw doorlopen. 

Criteria voor de selectie is een samenspel van een 
aantal factoren, zoals kosten, kwaliteit en beschik- 
baarheid. Het ligt daarbij voor de hand om de IP- 
verbinding als eerste keus te nemen. Maar als, bij- 
voorbeeld door overbelasting, de kwaliteit van de 
verbindingen over het IP-netwerk terugloopt, kan de 
PBX tijdelijk besluiten om nieuwe gesprekken over 
een alternatieve route te leiden. Bijvoorbeeld over 
een huurlijn of over het geschakelde openbare telefo- 
nienetwerk. 


Telewerkoplossing 

Telewerken refereert aan een breed scala van activi- 
teiten die worden uitgevoerd door een werknemer op 
afstand van het kantoor. Telewerken omvat altijd het 


SOPHO PBX 


Klant of collega 
maakt een gesprek 
naarde telewerker 


gebruik van telecommunicatie voor het contact tus- 
sen de telewerker en zijn collegae. Bovendien zal er 
in veel gevallen ook sprake zijn van één of andere 
vorm van elektronische informatieverwerking. Het 
gevolg is dat de telewerker toegang moet hebben tot 
het telefonienetwerk van het bedrijf en vaak ook tot 
het datanetwerk, het intranet van het bedrijf. 


De configuratie voor de telewerkoplossing is weer- 
gegeven in figuur 9. Personen, bijvoorbeeld een klant 
of een collega, die een gesprek willen opzetten naar 
de telewerker, hoeven niet te weten waar de telewer- 
ker zich bevindt. In alle gevallen wordt de telewerker 
gebeld op zijn zakelijk telefoonnummer. De telewer- 
ker zelf heeft volledige controle over de locatie waar 
zijn gesprekken worden aangeboden. Dit kan zijn op 
zijn werkplek op kantoor (A), op zijn normale tele- 
foontoestel thuis (B), op zijn PC thuis (C) of elke 
willekeurige andere locatie, bijvoorbeeld in de auto 
onderweg naar een klant (D). De voordelen die IP- 
Telefonie inbrengt in deze oplossing liggen in de 
configuratie bij de telewerker thuis. Vaak wordt ge- 
bruik gemaakt van een modemverbinding om de PC 
te koppelen met het bedrijfsdatanetwerk. Op dat 
moment is de thuislocatie telefonisch onbereikbaar 
voor normale gesprekken. Maar niet voor IP- 
Telefonie gesprekken, welke immers parallel met 
andere applicaties over de bestaande modemverbin- 
ding kunnen worden gevoerd. Tevens betekent dit 
een reductie van de gesprekskosten, omdat de tele- 
werker voor zijn uitgaande gesprekken geen extra 
lijn nodig heeft, maar gebruik kan maken van de toch 
al openstaande dataverbinding. 

De op IP-Telefonie gebaseerde telewerkoplossing 
bestaat uit een IP-Telefonie Gateway gekoppeld aan 


(0) 
Andere lokatie, 


@ 
Werkplek thuis 
PC-telefoon 

Á 


(@) 
Werkplek thuis 
normale telefoon 


Figuur 9: Telewerkoplossing over een IP-netwerk 
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of geïntegreerd in de bedrijfs-PBX voor de verbin- 
ding met het bedrijfsintranet. Via de IP-Telefonie 
Gateway zijn telefoongesprekken mogelijk tussen de 
PBX en de PC van de thuiswerker. De telewerkappli- 
catie bestaat verder uit een Telewerker Server op be- 
drijfslocatie en een Telewerker Client op de PC van 
de telewerker. 

De Telewerker Server bevat een database waarin 
voor elke telewerker staat geregistreerd waar de ge- 
sprekken moeten worden afgeleverd. De telewerker 
kan zelf de inhoud van deze database veranderen 
vanaf zijn thuis-PC. Dat betekent dat de telewerker 
zelf kan bepalen waar hij zijn gesprekken ontvangt. 


De telewerkapplicatie bevordert de bereikbaarheid 
van de telewerker door te monitoren of de telewerker 
bereikbaar is op zijn thuis-PC. Met andere woorden, 
de Telewerker Server is op de hoogte wanneer com- 
municatie met de Telewerker Client mogelijk is. 
Daarbij registreert de Telewerker Server tevens het 
huidige IP-adres waarop de thuiswerker kan worden 
bereikt. Dat is nodig omdat de toekenning van IP- 
adressen gewoonlijk dynamisch is en per sessie ver- 
schillend kan zijn. 


Web Call Centers 

Contacten tussen klanten en bedrijven vinden steeds 
meer plaats via het Internet. Elektronische commer- 
cie (“e-commerce”) is een geweldige groeimarkt. 
Bedrijven publiceren hun productcatalogus met tekst, 
plaatjes en soms geluid, op het World Wide Web. 
Soms echter, zal een potentiële klant behoefte heb- 
ben aan meer informatie en wil dan spreken met een 
vertegenwoordiger van het aanbiedende bedrijf. Dit 
kan heel goed worden opgelost met de technologie 


Multimedia klanten 
op het World Wide Web IP-Telefonie 


Gateway 


WWW Server 


van IP-Telefonie. Met een simpele klik van de PC- 
muis op een veld in de Web-pagina kan een spraak- 
verbinding worden opgezet tussen de PC van de po- 
tentiële klant en het toestel van een verkoper. Deze 
oplossing biedt voordelen voor zowel de koper als de 
verkoper. Voor de koper is het vooral gemak. Zon- 
der de Internet verbinding te moeten verbreken, kan 
gesproken worden met de verkoper, terwijl de Web- 
pagina nog op het beeldscherm staat. Voor de verko- 
per betekent het een kostenbesparing, doordat de 
klant niet later nog een keer moet worden terugge- 
beld. Bovendien neemt de kans op verkopen toe, 
doordat de koper meteen reageert en niet zijn interes- 
se verliest door een gedwongen bedenktijd. 


In figuur 10 is de applicatie geïntegreerd met de 
Automatic Call Distribution (ACD) functionaliteit 
van de SOPHO PBX. De verkooporganisatie heeft op 
een WWW-server een Web-pagina, die geschikt is 
gemaakt voor de IP-Telefonie. Dat houdt in dat de 
Web-pagina een software module bevat, zichtbaar 
voor de gebruiker (de koper) als icoon op de Web- 
pagina. Wanneer de gebruiker op deze icoon klikt, 
zal er een IP-Telefonieverbinding worden opge- 
bouwd tussen de PC van de gebruiker en de PBX van 
de verkooporganisatie, via de IP-Telefonie Gateway. 
Om dit te kunnen doen, moet de PC van de klant zijn 
uitgerust met een applicatie voor IP-Telefonie. Deze 
applicatie wordt automatisch op de PC van de ge- 
bruiker geïnstalleerd zodra deze de betreffende Web- 
pagina inleest. 

De gesprekken die op deze manier de PBX binnen- 
komen, worden door de PBX gedistribueerd naar één 
van de beschikbare verkopers, de ACD-agenten. 
Eventueel kunnen oproepen in een wachtrij worden 
geplaatst als er momenteel geen verkopers vrij zijn. 


ACD agenten metnor- 
male telefoons en een 
PC in communicatie met 
klanten op het WWW 


Figuur 10: Call Center in verbinding met het World Wide Web (WWW) 
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Het telefoongesprek tussen koper en verkoper bij de 
Web Call Center oplossing is direct gerelateerd met 
de Web-sessie. Zodra de spraakverbinding tot stand 
is gekomen, zal de verkoper op zijn PC dezelfde 
Web-pagina te zien krijgen als de klant. Alle infor- 
matie die de klant invult op een elektronisch formu- 
lier in de Web-pagina (bijvoorbeeld een klanten- 
nummer) zal meteen te zien zijn op het scherm van 
de verkoper. Daarnaast kan de verkoper ook bepaal- 
de Web-pagina's naar de PC van de klant sturen, bij- 
voorbeeld om de klant op een ander product te wij- 
zen. 


De migratie van de PBX 


Gedreven door ontwikkelingen als Computer Telep- 
hony Integration (CTI) en IP-Telefonie, zal de inte- 
gratie van de bedrijfstelefonie- en bedrijfdatanetwer- 
ken steeds belangrijker worden. De PBX zal daar een 
belangrijke rol in spelen. Er is geen enkel apparaat 
dat zo is toegespitst op de telecommunicatiebehoef- 
ten van de zakelijke gebruiker als de PBX. Het aantal 


mt 
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Telefoonboek 
diensten 


Desktop productivity 


Telefoniste 


Figuur 11: PBX als TCP/IP router voor alle media 


ten (zie figuur 11). Daarnaast zal er een duidelijker 
scheiding komen tussen de verschillende rollen die 
een PBX vervult. Er zal een verdergaande ontkoppe- 
ling plaatsvinden van de basistransmissie-, schakel- 
en interfacefuncties enerzijds en de faciliteiten an- 
derzijds. Dit is geïllustreerd in figuur 12. Daarbij 
zullen de faciliteiten, geboden door de Telefonie Ser- 
ver, in hoge mate onafhankelijk zijn van de onderlig- 


IP-telefoons 


Figuur 12:De PBX migreert naar ontkoppeling van Schakelmatrix en Telefonie Server delen 


faciliteiten dat wordt geboden is aanzienlijk, evenals 
de aansluitmogelijkheden op andere soorten netwer- 
ken. Natuurlijk zal de PBX moeten aanpassen aan de 
nieuwe eisen. Meer en meer zal het TCP/IP protocol 
zijn intrede doen. Niet alleen voor zaken als IP- 
Telefonie, maar ook voor netwerkbeheer, CTI- 
applicaties, telefoonboekdiensten en andere facilitei- 


gende infrastructuur. Gebruikers van IP-telefoons 
kunnen daardoor over dezelfde faciliteiten beschik- 
ken als gebruikers van “normale” telefoontoestellen. 
Tevens zal de Schakelmatrix de toegang blijven vor- 
men tot een veelvoud van andersoortige netwerken. 
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Conclusies en slotopmerkingen 


IP-Telefonie in de bedrijfsomgeving brengt een aan- 
tal voordelen. Het meest zichtbare voordeel op de 
korte termijn is een reductie van de kosten op com- 
municatie en infrastructuur. In de nabije toekomst 
zullen de toegevoegde functies steeds belangrijker 
worden en het is te verwachten dat IP-Telefonie een 
grote invloed zal hebben op de verdere ontwikkeling 
van spraakcommunicatie. 


Op dit moment kan IP-Telefonie nog niet dezelfde 
kwaliteit spraak bieden als de traditionele telefonie- 
dienst, de Plain Old Telephony Service (POTS). 
Maar de ontwikkelingen gaan snel, er komen nog 
betere compressietechnieken en hinderlijke effecten 
als spraakvertragingen zullen steeds meer worden 
teruggedrongen. Daarnaast komen middelen be- 
schikbaar om de capaciteit van IP-netwerken te be- 
heersen. In intranetten, toegepast in zakelijke omge- 
vingen, is dat vaak nu al mogelijk. 


Gebruikers van IP-Telefonie willen kunnen beschik- 
ken over dezelfde faciliteiten als de gebruikers van 
de traditionele telefoondienst. Daarom zal de PBX, in 
een nieuw jasje, een grote rol gaan spelen in het be- 
drijfscommunicatiesysteem van morgen. 


Eise afkortingen 


Automatic Call Distribution 
Asynchronous Transfer Mode 

Computer Telephony Integration 

Digital Enhanced Cordless Telecommunications. 
File Transfer Protocol 

Hypertext Transfer Protocol 

Internet Protocol 

Integrated Services Digital Network 
International Telecommunications Union (Tele- 
‘communication standardization sector) 
Local Area Network 

Private Branch eXchange 

Plain Old Telephony Service 

Public Switched Telephony Network 
Registration, Admission, Status 
Real-time Transmission Control Protocol 
Real-time Transmission Protocol 

Simple Network Management Protocol 
Transmission Control Protocol 

Trivial File Transfer Protocol 

User Datagram Protocol 

World Wide Web 

Wide Area Network 
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Impressies van de 
Nerg Werkvergadering van 3 november 1998 over het onderwerp: 


Humanitair Ontmijnen. 


De deelnemers zullen deze bijeenkomst zonder 
twijfel als heel bijzonder hebben ervaren en ze zullen 
er nog vaak over hebben nagedacht. Gaandeweg de 
middag werd het namelijk voor de aanwezigen in het 
TNO Fysisch- en Elektronisch Laboratorium te Den 
Haag duidelijk, welk een gruwelijk probleem de 
aanwezigheid van landmijnen in vele delen van de 
wereld vormt. 

Na de opening door de waarnemend voorzitter van 
het NERG, ir. W. van der Bijl, werd de reeks 
voordrachten geopend door dr.ir. BJ. Kooij van de 
TU Delft met het onderwerp “Inversion of a buried 
object” (zie ook het artikel op pag.) . De spreker 
dompelde het gehoor meteen in het diepe door de 
theoretische achtergrond van de mijndetectie middels 
elektromagnetische velden te behandelen. 

Voor de algemene probleemstelling kan men uitgaan 
van de onderstaande schets: 


Zenders / ontvangers 


Medium 1 


Medium 2 Verborgen 
object 


Het medium 2 stelt de aardkorst voor, waarin een 
inhomogeniteit voorkomt in de vorm van een 
onder de oppervlakte begraven voorwerp. De heer 
Kooij behandelde een methode om de 
materiaaleigenschappen van het begraven 
voorwerp te reconstrueren op basis van metingen 
aan de gereflecteerde signalen. 

Als ‘rode draad’ in het betoog kon ik het volgende 
optekenen. 

Er wordt uitgegaan van het optredende contrast in 
de elektrische veldintegraal van het object en z’n 
omgeving. Gemeten wordt in feite de per- 
mettiviteit van het object ten opzichte van de 
permettiviteit van de omgeving. Daartoe wordt het 
object aangestraald met een serie opeenvolgende 
excitaties waarvan de eigenschappen volledig 
bekend zijn. 

Optredende polarisatie-effecten rondom het object 


kunnen in de gereflecteerde signalen worden 
vastgesteld. De gradienten worden in verschillende 
richtingen vastgesteld. De patronen die daarbij 
ontstaan zijn afhankelijk van de vorm van het object 
. cilindrisch c.q. rechthoekig). De gemeten 
gradient functies kunnen nu worden vergeleken met 
berekende waarden. 

Vervolgens toonde de heer Kooij enige resultaten 
van uitgevoerde experimenten en beantwoordde hij 
nog een aantal vragen. 


De voordracht van prof.dr. habil. H. Rothe, 
verbonden aan de Universiteit van Hamburg, was 
getiteld: “Advanced Pattern Recognition”. 
Patroonherkenning vormt een belangrijk onderdeel 
van de moderne wijze van oorlogvoering. De 
situatie zoals in de tweede wereldoorlog, waar bij 
bombardementen veel bommen hun doel misten, 
ligt reeds ver achter ons, Maar in die periode ligt 
wel het beginpunt van de ontwikkeling van beeld- 
en signaalverwerking in het militaire toepassings- 
gebied. Prof. Rothe behandelde in kort bestek de 
methoden en technieken zoals die heden ten dage 
bij machinale patroonherkenning toegepast worden. 
Hij maakte duidelijk dat de keuze van geschikte 
kenmerken van het te herkennen object van 
doorslaggevende betekenis is. 

Voor de detectie van landmijnen werd door prof. 
Rothe een methode aangegeven, die hij omschreef 
als: ‘laser backscattering of the surface’. De 
onderzochte aardoppervlakte wordt bij deze 
methode aangestraald met 500 optische fibers. 

De signaalverwerking vindt plaats met een 
kunstmatig neuraal netwerk. Voor het “inleren’ (en 
dus ook tijdens het classificatieproces) worden 5 
kenmerken gebruikt. De c’ catie vindt plaats op 
basis van 4 verschillende algorithmen. Er zijn 
experimenten uitgevoerd op de _testlocatie 
Ft.Belvoir (U.S) en daarbij werd elke 1.5 inch van 
de locatie aangestraald. De einduitkomst van elke 
meting werd binair aangegeven, dus: wèl of geen 
landmijn aanwezig. Van alle positieve uitkomsten 
( landmijn aanwezig) was 90 % correct , dus 10% 
foutief. 

Maar, zoals later op de middag nog zou blijken, is 
met die cijfers nog niet alles gezegd over de 
praktische bruikbaarheid van een systeem. 
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Met name gaat het ook om de vraag, in hoeveel 
gevallen een landmijn ten onrechte niet wordt 
waargenomen. Bij de uitgevoerde experimenten op 
de testlocatie was dat < 1 %. 


De voordracht van ir. R.V. de Jongh ( TU Delft) had 
als titel: “The application of ground-penetrating radar 
(GPR) in the detection of buried landmines’. De heer 
De Jongh is verbonden aan het International 
Research Centre for _Telecommunications- 
transmission and Radar (IRCTR). 

Er bevindt zich een groot aantal landmijnen onder 
het aardoppervlak, naar schatting meer dan 100 
miljoen stuks. Dat vormt een wereldwijd probleem, 
omdat het om ca. 60 verschillende landen gaat. 

We spreken hier over humanitair ontmijnen; dat in 
tegenstelling tot militaire acties voor de eliminatie 
van mijnen. 
Globale eisen, te 
ontmijnen: 

— _ betrouwbaarheid van de detectie 99,9 % 

-__een redelijke snelheid bij terreinonderzoek. 
Conventionele systemen voor landmijn detectie, 
zoals prikstok, metaaldetector, speurhonden etc. 
voldoen hier niet aan. Daarom wordt er nu gewerkt 
aan de ontwikkeling van meer geavanceerde 
systemen, waaronder GPR. Het antenne-ontwerp 
vormt daarbij een speciaal probleem. Verdere 
bijzonderheden daarover vindt u in het artikel van ir. 
de Jong op blz. van dit nummer. 


stellen aan systemen voor 


Prof.dr.ir. Zwamborn gaf in zijn voordracht een 
nadere toelichting van de benadering die door TNO- 
FEL werd gekozen. Men gebruikt ultra brede band 
GPR voor de detectie en classificatie van objecten 
die zich onder het aardoppervlak bevinden. 
Classificatie is veel lastiger dan detectie. Ten 
behoeve van de technische uitvoering van de 
experimenten heeft men enige daarvoor belangrijke 
componenten gekocht: 
-_een generator voor het opwekken van zeer korte 
pulsen; 
- een Impulse Radiation Antenne voor het 
frequentiegebied van 200 MHz — 3 GHz; 
de diameter van die antenne is | meter. 
Bij het thans lopende project richt de aandacht zich 
op de ontvangst en synthese van het gereflecteerde 
veld. Met name wordt gewerkt aan ‘preprocessing’ 
en “feature extraction’ . 


Namens het Ministerie van Defensie werd het project 
‘“Humanitair Ontmijnen 2000’ toegelicht door ing. J. 
Roozenboom. Het project werd in 1997 gestart en 
heeft tot doel het ontwikkelen van apparatuur en 
methoden voor de detectie en verwijdering van 


landmijnen. Het past in internationale activiteiten, 
die onder auspicien van de Verenigde Naties 
aangestuurd worden. 

De door de V.N. gehanteerde eisen zijn: 

-__ 99,6 % betrouwbaarheid van de detectie; 
10.000 m° per dag te onderzoeken; 

- eenvoudige bediening. 

Men onderzoekt ook de eigenschappen van 
alternatieve technologieën, zoals microbiologische 
systemen, bio-sensoren, systemen waarbij neutronen- 
activiteit bepaald wordt etc. en ook GPR. 

Daar elk van die systemen z'n karakteristieke 
eigenschappen heeft, wordt ook onderzocht welke 
voordelen het combineren van verschillende 
systemen heeft. 

Van het verdere onderzoek aan GPR systemen 
verwacht men: 

—_optimalisatie van de bestaande technologie; 

- een lage kans op vals alarm; 

- ‘stand-alone’ gebruik; 

- identificatie van objecten. 


Als afsluiting van de bijeenkomst werd er een open 
discussie gevoerd, waaraan met name door de 
aanwezigheid van mevrouw Sylvia Martens van Pax 
Christi een bijzonder accent werd gegeven. 
Mevrouw Martens heeft enige jaren in Angola 
gewerkt, waar het probleem van de landmijnen al 
vele slachtoffers heeft gemaakt. Om goed voorbereid 
naar dat gebied te gaan, had mevrouw Martens 
vooraf een opleiding bij Defensie gevolgd. 

Voor de toehoorders was hetgeen mevrouw Martens 
vertelde bepaald niet opwekkend; het was “back to 
earth”. 

De harde praktijk is, dat in een land als Angola onder 
het aardoppervlak zeer veel levensgevaarlijke 
objecten van zeer verschillende aard zijn begraven: 
granaten, mortieren, landmijnen etc. Bovendien 
liggen er veel z.g. trip-lines. De gevaarlijke objecten 
liggen op zeer verschillende diepten, soms wel tot 9 
meter diep. Volgens de ervaring van mevrouw 
Martens zijn er veel “useless machines’die men inzet 
voor detectie. Het is namelijk heel wat anders om in 
die omgeving een gebied ‘gezuiverd’ te verklaren 
dan een goede score te behalen op een testlocatie 
waar zich alleen ‘echte’ landmijnen bevinden. Ook 
werd duidelijk, van welke betrekkelijke waarde de 
herkenning van een object is. In Angola zijn Coca- 
Cola blikjes gebruikt om landmijnen van te maken. 
In die omstandigheden hebben de bestaande 
definities over betrouwbaarheid weinig waarde. 

“Go to Angola and test your machines”. 


Arie Spanjersberg 
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NONLINEAR INVERSION OF A BURIED OBJECT 


B.J. Kooij 
Laboratory of Electromagnetic Research, 
Centre for Technical Geoscience, 
Delft University of Technology, 
PO. Box 5031, 2600 GA Delft, The Netherlands. 


Abstract 


A method for reconstructing the material properties of a buried bounded 
inhomogeneous object from measured scattered field data at the surface of the Earth's 
interface is presented. This work extends the method previously developed for the 
homogeneous TE-case to the more complicated case of a buried object. In the TE-case, 
the magnetic field is polarized along the axis of an inhomogeneous cylinder of 
arbitrary cross-section and the corresponding integral equation contains derivatives of 
both the background Green's function and the field. The non-linear inversion based 
upon the modified gradient method can be formulated as an electric field integral 
equation for the two transversal components of the electric field. The integrand is a 
product of the background Green's function, the contrast and the electric field vector, 
however in the case of a buried object the background Green's function is the one 
pertaining to a two-media configuration. The derivatives are operative outside the 
integral. In this paper the latter formulation will be taken as point of departure to de- 
velop a nonlinear inversion scheme. In the inversion scheme the positivity of the 
material parameters has to be ensured in order to obtain a convergent inversion 
scheme. For the inversion scheme the contrast source inversion scheme has been used. 
Numerical results are presented in which the reconstruction of the contrast is shown 
and compared in the case of TE excitation as well as TM excitation. 


1 Introduction 


We present a method for reconstructing the material properties of a bounded inhomogeneous 
object buried inside a half-space from measured scattered field data at the surface of the half 
space, see Kooij and van den Berg (1). The present paper extends the homogeneous TE-case by 
to the more complicated case of an object buried inside a half space. In the TE-case, the magnetic 
field is polarized along the axis of an inhomogeneous cylinder of arbitrary cross-section and the 
corresponding integral equation contains derivatives of both the background Green's function for 
the half-space and the field. This problem was investigated by Lambert et al. (2) using the scalar 
magnetic field integral equation. However, the integral equation can also be formulated as an 
electric field integral equation for the two transversal components of the electric field, see Kooij 
and Lambert (3). The derivatives are operative outside the integral operator. These derivatives 
can efficiently be integrated by using rooftop functions in the discretization procedure, see 
Zwamborn and van den Berg (4). The method is based on the modified gradient technique which 
is used by several authors in the same context, i.e. Kiemman and van den Berg (5),(6),(7), Lixin et 
al. (8). In this paper the electric field integral equation will be taken as point of departure to 
develop a contrast source inversion scheme using the modified gradient method. 

2 Integral Representations 


Consider the scattering object to be an inhomogeneous lossy dielectric cylinder of arbitrary cross- 
section D imbedded in a half-space with permittivity eh. The incident field consists of 
electromagnetic waves emitted inside a vacuum domain above the interface at x9 = 0 of the half- 


space with the magnetic vector polarized along the cylinder axis. To reconstruct the complex 
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permittivity of an unknown object from a knowledge of the scattered field we assume that the 
object is illuminated successively by a number of different excitations. For each excitation, time 
factor is exp(-i@t), we assume that the scattered field is measured exterior to the scatterer above 
the interface of the half-space. The complex permittivity of the buried inhomogeneous dielectric 
object is given by 

E(X) =E Ert iO/0, 

1) 


where €‚ denotes the relative permittivity of the object with respect to the permittivity of vacuum 
and 0 denotes the electric conductivity of the object. The incident electric field vector that is 
generated by a unit line source located at x= x* is denoted as E *“ =(E/", E‚“) and the incident 
magnetic field vector H“ has only the component H3”“. When we assume that that the incident 
field is the field that is present in the two-media configuration without the presence of the 


ine 


scattering object, the electric field vector E "“ =(E "5, E‚"°) inside the half-space is obtained as 


TYn 2e, 
a An (EoYn+EnYo 


x wol af nsi) +iYo 


in which yon = (Ron -£ 2, ko,h= @( Eo.nHo, } * are the wave numbers of the upper 
vacuum half-space, and the lower half-space, respectively. The electric field integral equation 
inside the half-space with the point of observation in D is obtained as 


en Je 


s 
x 
2 


AE @) 


E(&) - (la. + VV) AD =E""(@), xe D, 5) 


in which the vector potential A is defined as 


A= | Zag U fez) GE ere! eol) zee] (4) 


veD 


with 


Gr -x[ealtxal)= i jk EoYn—EnYo 


2 S= 2yileoyn +EnYo ) 
xexp Lia (x, =X,)+iY, (ll tal Jaa 5) 


and the normalized material contrast 


8)-E 
zi) © 
€, 
Once the scattered electric field vector, E““* = E - E“"“ is known in the contrasting domain D, 


ser 


the scattered magnetic field H ‚(we measure this field component) inside the vacuum half- 


space is obtained as 
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Vz. EB“ = H3"(x), D 


in which the vector operator V describes the transmission of the electromagnetic field into the 
vacuum half-space and where x” denotes the point of observation. The operator V is given by 


Vy.E“= Í G (ar a) (x'). Vond (2! Jav, ® 


xeD 


in which 


SEE zat S= 2E, 
ox -i] a Ì ár (ev, +e,v) 


xexplia (x1—x)+iy, ES +iy lesjlaa @) 


Note that the operator V is obtained from the unit line-source solution in the two-media 
configuration without the presence of the dielectric scatterer. 


3 Discretization Procedure 


We assume that the domain D is a rectangular domain with boundaries along the xj- and 
x-directions. We discretize the domain in a rectangular mesh. The mesh is uniformly 
spaced in the x,- and x,-direction. The rectangular subdomains Dy, m= 1... M‚ n= 1……, 
N with a width of Ax in the xj-direction and Axz have center points (2j:m, X2;:) where the 


contrast X is assumed to be constant with value yn Eq. (3) is discretized as 


Exam KiArimn B wma =Eloin Ke {1,2} (40) 


in which the function Bn, is the finite difference approximation of the gradient divergence 
operator that is acting on the vector potential Am, occurring in Eq. (3). The function Bimo in 
terms of the vector potential Aksm.n is then given by 


A a 


val pal gel 
where 


ad (AEP (E an DE ptn) an « 12) 


mara 


(2) (21) 
O mn.pa SC mn. 


1 5 
zr (4x Ax) é (s mai 0 an) 
(8 ons 8 one). 13) 


ol Ae) (BrranDO rat se 4) 


map gn 
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In Eqs. (10)-(14) Ekin = Ek ttms X2,a), Bomen = Bk (X1:mo,Xzn) and Agoma = AKQti;ms 2,0). 

Next we have to replace the continuous representation of the vector potential A by a discrete 
one. First we decompose the vector potential into a homogeneous part A* and reflected part A” 
according to 


A=A"+A" as) 


The Green function in the homogeneous part of the vector potential A° is non-regular in the 
domain of integration, while the Green function in the reflected part of the vector potential A’ is 
regular i in the domain of integration. In order to cope with the singularity of the Green function in 
AÌ, we take also the weak form by replacing it by the circular mean of the vector potential. We 
integrate Ag" over a circular domain with center at the point (;m, X2;n). The radius of these 
patches is taken to be 


LAx=L min (Ax‚‚Ax,). 
2 2 


The results are divided by the surface area id (tAx) de 
We then may write 


h 
Alm A (tims an) 
fo Gtm ttr) (oa )Er (ern )dr drs, 
with Ke {12} (16) 


where we have interchanged the order of integrations, such that 


zaate 4x mtr) de, 


5’ Jk, Ax}1$) k, (er +3)? ] 


k, 
n; arend, + 22)» 4 JAx, 
ä i 17 
Bla), Ee 
DET nk,Ax 


when x, =x, =0. 


In fact, G is the mean value of the Green function over a circular domain with center at (x1, X,). 
After this weakening procedure, we are now able to compute the integrals of Eq. (16) 
numerically. In view of the functional properties of Ek we approximate the integral of Eq. (16) 
using a midpoint rule in the xj- and x,-directions, We then arrive at 


Aims=Am an DDG! ant Kim |X 2 lenen (18) 
minst 
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Aln =Ax Ax, YYe (en T Krim 2 Kain T Kain )z mn E kome. 
and mnl 


(19) 


with Ke {1,2} and for m = 1 …. ‚M andn = 1 N. Note that Ana are discrete 
convolutions in xj and x3, while A'k;m» is a discrete convolution in x, and a discrete 


correlation in lx» |, which can be efficiently computed by FFT routines, see Lesselier and 
Duchene (10). In order to avoid numerical problems in the FFT routines due to the branchpoints 
of % and #, when we have zero conductivity in one of the media, we introduce a small 
eonductivity loss in the corresponding medium such that the branchpoint is no longer located on 
the real path of integration. 


4 Forward Scattering Problem 


When we substitute Eqs. (18) and (19) together with (15) in Eqs. (11) - (14) and use the result in 
Eq. (10), we obtain a linear system of equations for Ex; When the material contrast 4 is 
known. This system of equations can be solved by a conjugate gradient scheme. In such a scheme 
we need the representation of the operator, the definition of the norm and inner product and the 
representation of the adjoint operator. The operator expression LE is directly obtained from Eqs. 
(10)-(15) and Eqs. (18) and (19) as 


2 M N 
(LE) ns =E Kira ki A Kip DN a GA À vonnis (20) 


for p= 1…, M andq= L…, N, and where Ag.» directly follows from Eq. (15) and Eqs. (18) and 
(19). 


The norm on D is defined as 
2MN pe 
VED -EID Bers eN 
vl pel gel an 


where the overbar denotes complex conjugate. The adjoint operator is defined through the 
relation 


(HLB) (LRE) 5 (22) 


Substituting the expression of the operator LE in the left-hand side of Eq. (22) and interchanging 
the various summations, the adjoint operator is recognized as 


(L*r Vesoa ST Kipa AT =Ckipa 
(23) 


for p= |. , M and q = 1…, N, in which 


C gina Chipo +C kina 4) 
with 
2 MAIN=L 
C Kipa =Ax Ax, DINER Trip: 2 Hog Kaigr )o Kipa” 
p=0q0 
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(25) 


and 
MINA _R 
Chips =Ar Ax DDG (arti 4 
p'=0q"=0 (26) 
where 
ie 2 M N « ) 
Ona ker mima t Dede, El a rem en 


vel mel n=l 


The latter expressions may be written explicitly as 
—i 2 
O ina =k Hina + (4 % ) ear hija Hhipng ) 


l „ 
+7 % Ax) s e Zip TT pan Taipan FT oip gn ). (28) 


[) Bp =kj Vaipa (Ax) bene 24 +Fipgn ) 
1 4 
ried Ax) ï (- Bip-har Tp TT vaer Fisprgn ): (29) 


Since, according to Eqs. (25) and (26) p' runs from O to M + 1 and q' runs from O to N + 1 
we set in Eqs. (28) and (29) 


rk:p.4= O for p=-1,0, M+1, M+2, Vq, 
FK:p.y= O for q=-1,0, N+1, N+2, Vp. (30) 


Note that C'x‚pg are discrete convolutions in xj and xs, while C'xg is a discrete convolution in 
xj and a discrete correlation in bel, which can be computed efficiently by FFT routines as 
mentioned in the previous Section. 


With these definitions we are now able to apply a conjugate gradient iterative scheme to 
solve the equations 


(LE) 4 =E Kpn P=leeM, q=l.N GI) 


Once the normalized error, 


Îre-e"|p 2) 
le=[p 


is small enough, the approximate solution of E is substituted in Eq. (7) to arrive at the scattered 
magnetic field at the receiver points x”, withr= 1,2... R.. 


ERR = 


5 Contrast Source Inversion Method 


We now assume that the inhomogeneous object D is irradiated successively by a number (s = 1, 
2,…, S) of known incident fields. To show this explicitly we rewrite Eq. (3) in terms of 
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the contrast source WO) =y E © 


WO_y LW = BE", on D, (33) 
where 
Z.M 
(s) EE Kv (s) 
(Lw ee =k, Akira +22 PA (34) 


for p =l…, M and q = 1…, N. Equation (33) is denoted as the OBJECT EQUATION that holds in 
the object domain D. 


Further, in the inverse problem H3 is known at the measurement points x'. We assume that 


all the measurement points are located in the data domain S. Hence, Eq. (7) is rewritten in terms 
of the contrast source as, 


v.w”= 3x”), on S, (25). 
is denoted as the DATA EQUATION that holds in the data domain S. 


The inverse problem is that of finding the material contrast x from Eqs. (33) and (35). For 
the inversion we construct sequences of the contrast sources W‚*, the fields E‚® and the 
contrasts ; to minimize the cost functional 


F=D' 5 [a®-v.w Ol b 36) 
sel 


2 |-1 
D 


5 S s -l 
and p'=(5", | Hi Î Í ;) i where the subscripts D and S are included in the norm || || in L2 


+ e” @O_w Oz Lw © | 


in which D° 5 [y E*“O 


to indicate the domain of integration. Next, we construct sequences {W‚® } and (%4}' fork = 
1,2…, in the following manner. We define the data error to be 


pP=HP-v.wld, G7) 

and the state vector to be W=z,ED-wWO, (38) 
where 

EO zE" O4 LW (9) 


If we assume that W”‚, and 4, are known, we can update W‚“ by 
Ip 
WOW ta, wb) 40) 


(e) 


The update directions w‚ “) are chosen as the Polak-Ribiere conjugate gradient directions, namely 


w{) =0, 
we) =dwO4y “wil for kel, GI) 
where 
S 
Naw aw Paw P)D 42) 
w ssl 


yhowls 
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In Eqs. (41) and (42) the function Öw,” is the gradient or Frechet derivative of the cost 
functional with respect to W‚” evaluated at W‚,® and 4, and (…) p denotes the inner 
product on L2(D). The Frechet derivative was obtained as 


awP=-D*V*p- or [- elzen) (43) 


in which the * symbol denotes the adjoint of the operator and the overbar denotes complex 
conjugation. 

2 |l 

5) 


InEa.(43)  DP=l bi 
The constant o4 can now be obtained from the cost functional F as 


Iz EE aal 


EC: Xol a Mew Oise DP „Del owbozaate | 


5 + Dre! DL a Lw) 


ogv 


where (.,.)s denotes the inner product on L2(S). The function E‚ 
(39) and (40) as 


)} (44) 
Ake 


©) can now be found via Eqs. 


EO =EO +a,Lwb, (45) 
and we seek to minimize the cost functional 


‚=de — 


For this cost functional we have not proven yet that the error is reducing in the iterative 
procedure, however in all our numerical examples it was observed that the error was reducing 
monotonically. The form of the cost functional presented in Eq. (46) is very suitable for easy 
implementation of a priori information or constraints on  . 


(46) 


We write X as 
xt! (47) 
Then in the absence of a priori information we find that Fp is minimized by choosing 


ZRD. EO) 
k EP (48) 


and 


(49) 
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However, when a priori information is available we find that Fp is minimized by choosing 


(50) 


and 


- 61) 


if x “and x “are positive. 


6 Numerical Examples 


In our numerical examples it was assumed that the unknown scatterer was located entirely within 
a test square of known dimension although knowledge of the precise location within the test 
square was not assumed. The convolutions and correlations that occur in our numerical scheme as 
well as the inverse Fourier transformation of Eqs (2) and (5) were carried out using FFT 
techniques which speed up the computation time significantly. We only present numerical 
examples in which the positivity of the permittivity or conductivity is ensured, Without the 
insurance of the positivity of these material properties our iterative scheme converges to a 
reconstruction with negative non-physical material properties. Probably due to the lack of data 
(we have only data from a limited viewing angle) we cannot simultaneously reconstruct the 
conductivity and permittivity, see Habashy et al. (11), Kooij and van den Berg (1). All 
computations were carried out on a dual 300 MHz Pentium II Workstation running Windows NT 
4.0. The FORTRAN code was compiled using a Digital Visual FORTRAN 5.0 compiler and the 
computation time for one iteration in the inversion in the examples of figure l and 2 was 1.8 
minutes. 


6.1 Example 1 


We illustrate in Figure 1 the TE-inversion as well as the TM-inversion of a positive conductive 
contrast with a numerical example, in which the measurements and excitations were taken on the 
same spot at 37 different positions at a distance of 0.125 m above the interface of the two media. 
The horizontal distance between the sources or receivers is set to O.l m. The complex relative 
permittivity of the upper medium is set to e,= 1.0 + 10° and the relative permittivity of the 
lower medium is set to &,= 2.0 + 10°i, The sources emit electromagnetic radiation at a frequency 
of 300 MHz. The test domain D, which is a square of à by A, where À is the wavelength in the 
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medium below the interface, is discretized into 61 x 261 subsquares. The center of this square is 
located 0.75 m below the interface. The configuration to be reconstructed consists of a square 
cylinder located in the center of the test domain D and is discretized into 23 x 23 subsquares 
where each subsquare is of dimension 11.592 mm x 11.592 mm and has a relative complex 
permittivity &,= 2.0 + 2.0 i 


In the examples, we assume that 37 stations (J = 37) are located equally spaced above 
the interface S with each station serving successively as a line source emitting radiation at 
a frequency of 300 MHz and all stations acting as receivers. The contrast x is normalized 
to a real quantity in the depiction of our numerical results in this example and is given by 


(no) ele, 


pd 
el’ 


62) 
In Figure 1 the actual normalized contrast x (0) is depicted in (a), while the results after 32 and 
128 iterations for TE- as well as TM-polarization are depicted in (b), (d) and (c), (e), respectively. 


6.2 Example 2. 


We illustrate in Figure 2 the TE-inversion as well as the TM-inversion of a positive permittivity 
contrast with a numerical example, in which the configuration is identical to the one in the 
example la. However, the square cylinder located in the center of the test domain D which is 
discretized into 23 x 23 subsquares where each subsquare is of dimension 11.592 mm x 11.592 
mm now has a relative complex permittivity €, = 4.0 + 10'Sí. 


The contrast % is again normalized to a real quantity in the depiction of our numerical results 
in this example and is given by 


_el)—e, 63) 


In Figure 2 the actual normalized contrast ME) is depicted in (a), while the results after 32 


and 128 iterations for TE- as well as TM-polarization are depicted in (b), (d) and (c), (©), 
respectively. 
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(b) 


(d) 


Figure 1: The reconstruction of a conducting square cylinder with a bounded and 
real contrast Xx!” depicted in (a). The figures (b) and (d) represent a bounded 
reconstruction with TE-scattering after 32 and 128 iterations. respectively. The 


figures (c) and (e) represent a bounded reconstruction with TM-scattering after 32 
and 128 iterations, respectively. 
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Figure 2: The reconstruction of a square cylinder with a positive real permitivity 
contrast Xx“) depicted in (a). The figures (b) and (d) represent a bounded recon- 


struction with TE-scattering after 32 and 128 iterations, respectively. The figures 
(c) and (e) represent a bounded reconstruction with TM-scattering after 32 and 128 


iterations, respectively. 
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7 Conclusions 


We have developed a nonlinear inversion method for the two-dimensional problem of scattering 
of an electromagnetic wave with the magnetic field vector parallel to the cylindrical structure. 
The electric field vector has two components in the transversal plane (TE-case). Comparison of 
the numerical results calculated in the TM-case with those calculated in the TE-case shows that 
the results are different in these cases. It seems that inversion in the TE-case yields better results 
for penetrable objects while the TM-case seems to give better results for impenetrable objects. 
The main reason for these differences is that in the TM-inversion scheme the E3 component is 
used as the unknown quantity which is continuous at all edges (tangential incidence) of the 
contrast. While in the TE-inversion scheme both the Ej and Ez are the unknown quantities of 
which one has normal incidence and one has tangential incidence at an edge of the contrast. And 
since we have a limited view of the object, i.e. a limited set of measurement and excitation points 
present at one viewing side on the surface of the interface between the two media, the results for 
the TM-inversion and TE-inversion show large differences. 


The present solution of the TE-case is an intermediate step to the full threedimensional 
problem where all the three components of the electric field play their role. In the forward 
problem, the step from the two-dimensional TE-case to the full three-dimensional case is rather 
trivial. All fundamental difficulties of the three-dimensional case are the same in the TE-case. 
Only, in the three-dimensional case, there is extra book keeping in writing the computer code, 
and obviously for realistic configurations the computation time increases substantially, see 
Zwamborn et al. (12). It is expected that the same will be true for the inverse problem. Hence, it 
is a prerequisite to solve first the inverse problem in the TE-case. For the latter case, we have 
described the method based on the contrast source inversion method; this avoids the need of a full 
solution of the forward problem in each iteration of the inverse problem. With this feature it 
seems that a three-dimensional inverse problem can be handled with the present-day computer 
power. 
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Abstract 

One of the most critical hardware components for the performance of a ground 
penetrating radar (GPR) is the antenna system. Most GPR systems operate in the 
vicinity of the ground. In this case, the important parameters determining the 
system performance, such as antenna beamwidth, operational frequency and 
transmitted wavelet, depend on the properties of the ground. These properties 
change for different locations and weather conditions. The dielectric properties, 
roughness of the ground interface and stratification of the subsurface will 
influence the optimal antenna design. 

In a joint research project, IRCTR and RTI are developing new GPR antenna 
concepts. Aim of the developments are the improvement of GPR performance. 
The design of GPR antennas differ from the classical approach. Therefore, this 
paper is started with a review on the most important design principles for GPR 
antennas. To verify the antenna design, a GPR antenna measurement set-up is 
essential. Because little information on this topic could be found in literature, we 
developed our own measurement set-up and procedures. In this paper the first 
experimental results obtained from the set-up are shown. In the last part of this 
paper some GPR antenna types are discussed based on the experiments. Some 
new approaches in antenna design are presented as well. 


Keywords: GPR antenna design, GPR antenna measurements, short-pulse 
radiation, TEM horn, shielded broadband dipole. 


1. Introduction 


At Delft University of Technology, a large program is 
initiated to improve the current generation ground 
penetrating radar (GPR) systems. This project is 
financially supported by the National Technology 
Foundation (STW) and is named “Improved Ground 
Penetrating Radar Technology”. Three groups of the 
Delft University participate, the Centre of Technical 
Geoscience, the Electromagnetic Laboratory of the 
faculty of Information Technology and Systems, and 
the International Research Centre for 
Telecommunications-transmission and Radar 
(RCTR). 


The main task for IRCTR is the development of a 
transportable ground penetrating radar. The radar is 
based on transmission and reception of short-pulse 
time-domain signals. The antennas are located on or 
just above the ground. The work is done in 
cooperation with GeoZondas, Vilnius, Lithuenia and 
Academician A.L. Mints Radiotechnical Institute 
(RTI), Moscow, Russia. GeoZondas is responsible 
for the transmitter, receiver, controlling and 
processing software. RTI is working on the antenna 
design. 


The first part of this paper reviews the important 
design principles of the GPR antennas. The second 
part describes the GPR measurement range. The last 
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part of the article discusses our progress in the GPR 
antenna development. 


2. Overview design principles for GPR antennas 


The design of GPR antennas differs from the classical 
approach. It is therefore essential to formulate the 
relevant GPR antenna design principles. We restrict 
ourselves in the overview to antennas for short-pulse, 
time domain GPR, operating near the ground. 


2.1 Time domain antenna performance 


The desired antenna performance is to transmit and 
receive short duration, time domain waveforms (in the 
order of a few nanoseconds) with a little tail. The 
short duration of the pulse is used to obtain sufficient 
range resolution (Daniels). The tail of the antenna 
signal must be minimized to prevent masking of 
targets by nearby object reflections. The duration of 
the time domain antenna pulse is the trade-off 
between range resolution and penetration depth. 


To transmit short duration time pulses, an ultra- 
wideband antenna is needed. However, the bandwidth 
is not the only condition. A linear phase characteristic 
and constant phase center are also required. For 
example, spiral and log-periodic antennas have a high 
operational frequency range, but the time domain 
antenna response have a long duration. Note, that 
synthetic pulse reconstruction techniques are not 
attractive in time domain GPR systems due to the 
unknown propagation conditions in the ground. 


To discuss the time domain antenna performance we 
will introduce some terms and definitions. The class 
of antennas with ultra-wideband and linear phase 
properties is called transient antennas [Foster, 
1992]. The term transient is more closely related to 
time domain performance. The time domain antenna 
response is divided in two parts; the main pulse and 
the tail. In the main pulse we distinguish two regions 
(see figure 1). The transient region is resulting from 
the direct radiation of the excitation pulse. The 
resonance region is caused by reflections in the 
internal antenna structure (e.g. end of dipole, antenna 
shield). The tail of the antenna is caused by the non- 
ideal antenna characteristics. The tail is specified 
according some decaying law (e.g. exponential). 

The design goal for the time domain antenna 
performance is removal of the resonance region and 
minimizing the tail. Examples of transient antenna 
performance are the infinite dipole antenna and TEM 
horn. 


transient resonance 
region 


region 


main tail 
pulse pulse 


Figure 1: Definition of time domain antenna 
characteristics 


Major difficulty in the practical design is the removal 
of the resonance region. Due to the finite dimensions 
of the radiating element (e.g. a bow-tie), the pulse will 
reflect on the boundaries and result in multiple 
excitation points and times. A method of reducing this 
problem is resistive loading. In literature various 
solutions are proposed, ranging from resitive coating 
to the use of resistors on the boundary of the radiating 
element [references!]. 


2.2 Reduction ground influence 


For GPR systems located on or just above the 
ground, the current distribution on the antenna will be 
influenced by the ground. The radiation properties 
depend on ground conditions. The ground conditions 
change for different soil types, weather conditions 
(humidity) and surface roughness. Changing 
conditions restrict the interpretation of GPR data. 
Effects, like antenna ringing due to the antenna- 
ground interface, can make GPR measurements 
useless. 


The design goal for the GPR antenna is to reduce the 
influence of the ground on the antenna performance. 
A method to reduce the influence is to elevate the feed 
point of the antenna from the ground interface [ref]. 
Dielectric embedding of the feed point can lead to 
further improvements. The dielectric embedding 
reduces the contrast between the antenna and ground 
(thus reducing the influence) and increases the 
electrical length between feed point and ground 
interface. 
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2.3 Radiation characteristics 


In most GPR systems portability of the system is 
required. Because of the portability requirement and 
the relative low operating frequency range (most 
systems operate in the range 50 Mhz - 5.0 Ghz), GPR 
antennas have broad polar radiation patterns (low 
gain). Knowledge on exact antenna radiation is 
important for interpretation and processing of GPR 
data (eg. for migration). However, the classical 
characterisation of antennas by its far-field radiation 
pattern is of little use in GPR antennas because of 
two reasons. First, objects can be located in the near- 
or intermediate region of the antenna. For processing 
of this data, the far-field pattern cannot be used. 
Second, in most cases the ground will have a layered 
structure. Due to the layers, the actual far-field of the 
antenna is depending on electrical properties, 
structure and thickness of the layers. 


A more practical method of evaluating the radiation 
characteristics is to define the EM field on a planar 
interface just under the ground interface. Using near- 
to-far field transformation techniques in time domain 
[Yee etal. 1991), the radiation pattern for each 
specific ground can be determined (for example using 
EM modeling). In section three of this paper a 
measurement method is proposed to characterize 
GPR antenna radiation. 


2.4 Antenna shielding 

The desired GPR antenna radiation is directed 
towards the ground. Radiation in the upper halfspace 
(air) must be minimized to reduce the influence of 
interference signals (eg. tv, mobile phones) and 
remove undesired reflections from targets above the 
ground. 


To minimize the radiation to the upper halfspace 
antenna shielding can be used. Some antenna types 
are already shielded (like the horn antenna), while 
others need an additional design effort (like the dipole 
antenna). Major difficulty is that the shielding must 
not reduce the time domain performance. 


3. GPR antenna measurement set-up 

Antenna measurements are an essential part of the 
design procedure. However, little information can be 
found in the literature on the specific topic of GPR 
antenna measurements. Therefore, we are developing 
our own measurement set-up and procedures. 


3.1 Test site 


For GPR antennas and systems it is essential to 
measure the characteristics in a well known test site. 
Together with FEL-TNO, the Delft University of 
Technology constructed a GPR test range, located in 
The Hague. The dimensions of the range are 
10x10x3 m. The large dimensions reduce sidewall 
reflections which interfere with the measurements. 
The test range is filled with ‘homogeneous’ sand. To 
sustain controlled conditions, the test range is 
shielded from external influences. The test-range is 
separated from the surrounding ground by wooden 
walls, covered with a water proof plastic coating. 
Weather conditions are excluded by using a tent. 


3.2 Measurement set-up 


As was discussed in the previous section, the far-field 
pattern of the antenna is of little use. It is more 
practical to define a surface under the ground and use 
transfomation algorithms to determine the radiation in 
the ground. From the near-field theory (in air) we 
know that for the determination of the antenna 
radiation, it is sufficient to measure on a specified 
surface two independent field components (out of six 
components). Additional research is needed to 
determine additional constraint for measuring the field 
in the ground. 


To realize the measurement set-up, an EM-field probe 
is buried just below the ground interface. The probe 
measures the local time domain EM field. The 
antenna under test (AUT) is moved on a planar 
surface above the ground. In discrete positions, time 
domain waveforms are measured and stored. For the 
planar movement of the AUT an X-y-z scanner 
(1.5x1.5x0.3 m) is developed. 

Measurement set-ups for other GPR characteristics 
(such as ground influence) are not described in this 
paper. Additional research on this topics is still 
necessary. 
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Figure 2: GPR antenna radiation measurements. 
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3.3 Measurement equipment 


Different pulse generators can be selected to excite 
the AUT. In table 1 the main pulse characteristics are 
summarized. The time domain signals are measured 
with a sampling oscilloscope (K2-63, table 2). 


Table 1: Pulse generator characteristics. 


Waveform | Duration Amplitude 
[ns] [Vv] 

monocycle | 0.8 [0% level] +25 
‘Ins’ | | 
monocycle | 2.6 [0% level] | + 100 
‘2.5ns’ 
monoeycle 5.0 [0% level] +250 
‘Sns’ | | 

Limpulse 1 [0.2 [50% level} | 40 


impulse 2 0.4 [50% level 50 


Table 2. Characteristics sampling oscilloscope 


Parameter Value 

Number of channels | 4 

Bandwidth 2 x(1+6 GHz) 

2 x(1+18 GHz) 

Noise (RMS) < 1.0 mV (1+6 GHz) 
- 4.0 mV (1+18 GHz) 

Max. input voltage 

Accuracy 


A/D converter 
Number of samples 


256, 512, 1024, 2048 


Max. averaging 256 
Min.sampling time 0.1 ps 
Jitter <20ps 


Measurementrange | 10.0 ps=l0s 


3.4 EM field probe 

The probe used in the first experimental antenna 
measurement is a home-made H-field sensor (loop 
antenna with inner diameter 30 mm (Goedbloed, 
1991). This simple probe antenn is used to 
demonstrate the feasibility of the measurement set-up. 
The probe is buried at depth 25 cm. The set-up is 
tested with two different AUT; an identical loop 
antenna and a TEM horn designed for Ins monocycle 
transmission. Figures 3 and 4 show the measured 
responses. 


Probe response 


02 } } 


0.1 
Amplitude 


Time [ns] 


Figure 3: AUT is loop antenna. 


‘Ins monocycle TEM horn response 


Figure 4: AUT is TEM horn 


The main pulse duration and tail for the TEM horn 
antenna signal is similar to the response of two 
identical TEM horns in the air. Furthermore, the 
response of the two loop antennas show no large 
resonances. Additionally, the signal strength of the 
different antenna responses is sufficiently large. We 
concluded that the measurement set-up used is 
feasible. Further analysis of the probe characteristics 
is not done, because in the final measurement set-up 
other probe antennas will be used. 


4. Analysis of different antenna concepts 

Aim of the research is to develop improved GPR 
antenna concepts. In this paper the first experimental 
results are described. Two antenna types are 
discussed: the TEM horn and the dipole antenna. 
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4.1 TEM horn 

An antenna type which is known to be good for time 
domain pulse radiation is the TEM horn. An air filled 
TEM horn is developed for the radiation of the 1 ns 
monocycle pulses (see figure 5). The antenna signal is 
shown in figure 4. 
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Figure 5: TEM horn antenna for Ins monocycle 
excitation 


For the TEM horn response the main antenna pulse is 
short in time and the tail is small compared to the 
main pulse amplitude. Note that from figure 4 the 
absolute antenna characteristics cannot be 
determined, because no probe correction is aplied. 


The TEM horn have also disadvantages. First, the 
crosscoupling is large. Second, the TEM horn is not 
shielded. Currently we are investigating improve- 
ments to this antenna concept for use in a GPR 
system. The first step is dielectric filling of the TEM 
horn. The EM waves propagating will be more 
contained in the dielectric embedding. This may 
reduce the crosscoupling between the antennas. 
Furthermore, the dielectric filling reduce the contrast 
with the ground, resulting in a smaller ground 


influence. Finally, due to dielectric filling the 
dimensions of the antenna can be reduced. 


4.2 Shielded broadband dipole 


The second antenna type investigated on the GPR test 
range is the ‘shielded broadband dipole’. This concept 
resulted from a compromise between two GPR 
antenna requirements. First, a shielded antenna 
structure. Second, small antenna dimensions. The 
concept is the placement of a broadband dipole in a 
metal casing. The metal case act as a waveguide and 
the dipole is used for waveguide excitation. Because 
the propagation path in the waveguide is short, the 
dispersion (typical for a waveguide) is expected to be 
small. The dipole arms have an elliptical shape 
(optimized for maximum flat frequency response). In 
the first prototypes, the dipole was placed in air. The 
first antennas showed large antenna ringing (see 
figure 6). 
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Figure 6: First prototype shielded broadband dipole. 


The time domain antenna response in figure 6 is not 
accepatable for the GPR system. By covering the 
sidewalls with absorbing material, the antenna ringing 
is reduced. The next step is the dielectric embedding 
of the antenna to reduce the contrast with the ground 
(see figure 7). The response of the improved antenna 
is given in figure 8. 


Figure 8 showed a large improvement comparing to 
figure 6, especially for the tail. However, further 
improvements can be made in the main pulse. In the 
near future experiments will be done with resistive 
loading of the dipole and better absorbers on the 
dielectric embedding. Another idea which is currently 
investigated is the use of lossy dielectrics 
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dielectric, 
dipole arm 


Figure 7: Top view dielctric embedded broadband 
dipole. 


in the antenna. The feed point is embedded in the 
dielectric. The dipole radiates to all directions. For 
the antenna performance it is better that all energy 
outside the direct path to the ground is absorbed. This 
can be done by inserting losses in the medium (carbon 
or ferrite particles). Because the electric permitivity 
of the medium is constant, no reflections are 
introduced. 
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Figure 8: Measured response embedded broadband 
dipole 


Conclusions 


In this paper the design of GPR antennas is described. 
First, a review of the main design principles is given, 
because GPR antenna design differs from the 
classical approach. Most important principles are a 
good time domain performance (removal of 
resonances and minimize the tail) and minimal ground 
influence. Second, a practical GPR measurement set- 
up is described. Measurements are essential in the 
antenna design. In the GPR test site we measure the 
radiated time domain field on a planar surface in the 
ground. The time domain antenna performance can be 
investigated for different excitation pulses. Third, two 
antenna types are discussed; the TEM horn and the 
‘shielded broadband dipole’. The TEM horn have 
good radiation characteristics, but large 
erosscoupling and is sensitive for interference signals. 
Additional work on the TEM horn is discussed. The 
‘shielded broadband dipole’ is designed to have a 
closed structure with small antenna dimensions. 
Intermediate results of the first and the latest 
prototype are shown. However, further improvements 
of this antenna type is expected. 


Voordracht _ gehouden _ tijdens de _467e 
werkvergadering. 


De auteur: 


31 Tijdschrift van het Nederlands Elektronica en Radiogenootschap deel 64-nr.1-1999 


René de Jongh is afgestudeerd in 1994 aan de 
Technische Universiteit Delft, Faculteit 
Elektrotechniek. Van 1994 tot 1996 is hij TWAIO 
aan dezelfde faculteit, bij het onderzoeksinstituut 
IRCTR. Onderwerp van de TWAIO studie is de 
ontwikkeling van een nieuwe meettechniek voor hoog 
frequente antennes. In deze methode worden antenne 
eigenschappen bepaald met behulp van korte 
tijdpulsen. Deze ontwikkeling wordt gedaan samen 
met het Russische Radiotechnical Institute (RTI). Ten 
behoeve van dit onderzoek wordt twee maanden 
doorgebracht op dit instituut in Moskou in de zomer 
van 1995, Na het afronden van de TWAIO studie in 
1996 wordt René de Jongh toegevoegd onderzoeker 
binnen het STW project “Improved Ground 
Penetrating Radar Technology”. Binnen dit project is 
hij verantwoordelijk voor de ontwikkeling van een 
grondradar systeem. Voornaamste onderdeel van dit 
systeem zijn de grondradar antennes. Deze antennes 
moeten in staat zijn om korte tijdpulsen te zenden en 
ontvangen en optimaal te werken in de nabijheid van 
de grond. Vooral deze laatste eis is moeilijk te 
realiseren omdat de grond eigenschappen sterk 
afhangen van de samenstelling, ruwheid van het 
oppervlak en de weersomstandigheden. 


Voortvloeiend uit zijn activiteiten bij het grondradar 
onderzoek is hij betrokken bij het specifieke 
onderzoek naar landmijn detectie met behulp van 
grondradar. 

De activiteiten met betrekking tot de korte puls 
meettechnieken en grondradar zijn verricht bij het 
International Research Centre for Telecommuni- 
cations-transmission and Radar (IRCTR). Dit 
instituut is onderdeel van de faculteit Informatie 
Technologie en Systemen (ITS) van de TU Delft. Eén 
van de onderzoeks gebieden van dit instituut is 
grondradar. IRCTR is verantwoordelijk voor de 
technologie ontwikkeling in bet STW project 
“Improved Ground Penetrating Radar Technology". 
Dit project, dat gestart is in 1996 en doorloopt tot 
2000, richt zich op grondradar toepassingen in het 
algemeen. Daarnaast wordt er voor het Ministerie van 
Defensie gericht onderzoek gedaan naar landmijn 
detectie met behulp van grondradar. Vanaf januari 
1999 begint IRCTR aan een groot onderzoeks project 
naar landmijn detectie met behulp van grondradar. 


Vanaf 1 oktober 1998 is René de Jongh werkzaam als 
technologisch ontwerper bij Hollandse Signaal 
apparaten B.V, bij de vestiging Delft. 


Tijdschrift van het Nederlands Elektronica- en Radiogenootschap deel 64-nr.1-1999 32, 


DESIGN OF AN ULTRA-WIDEBAND GROUND-PENETRATING RADAR SYSTEM USING 
IMPULSE RADIATING ANTENNAS 


J.B. Rhebergen and A.P.M. Zwamborn D.V. Giri 

TNO Physics and Electronics Laboratory, Pro-Tech, 

The Netherlands. United States of America. 
ABSTRACT 


At TNO-FEL, one of the research programs is to explore the use of ultra-wideband 
(UWB) electromagnetic fields in a bi-static ground-penetrating radar (GPR) system 
for the detection, location and identification of buried items of unexploded ordnance 
(e.g. land mines). In the present paper we describe the current status of the 
development of this system. 


The UWB ground-penetrating radar system is designed to operate in the frequency 
band from 200 MHz to 3 GHz and uses impulse radiating antennas (IRAs) as 
transponders to radiate and receive very short electromagnetic pulses from a short 
distance above the soil. The receiving IRA is similar to the transmitting IRA and is 
connected to a receiving unit which comprises of an attenuator, time gate switch, 
trigger delay generator and a sampling oscilloscope (see figure 1). A personal 
computer is used to control the equipment. In order to perform controlled radar 
experiments, a full-size experimenting and testing facility has been erected on the 
premises of TNO-FEL. 


SYSTEM DESIGN the main interest 

is in objects varying in size from about 5 cm up to 
ca. 75 cm. In air this corresponds to wavelengths of 
1.5 m to 0.1 m. Consequently the frequency range of 


interest extends from 200 MHz up to 3GHz. 


To design an UWB GPR system for detection, 
location and identification of land-mines the 
following points have to be successfully combined. 


« generation of UWB EM field (i.e. source) 

e radiation of UWB EM field (i.e. antenna) 

e __characterization of the soil and air/earth 
interface 

e reception and synthesis of scattered fields 

e__ signal processing techniques and identification 
of objects 

e _ controlled experiment/test environment. 
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This paper concentrates mainly on the system 
design. The signal processing aspects are, for the 
time being, saved for a future follow-up paper 


although some preliminary work has taken place [1]. 
An overview of the overall system concept is shown 
in figure 1. The target volume is illuminated by an 
ultra-wideband electromagnetic field, with an angle 
of incidence chosen to maximize energy 
transmission into the ground, whilst minimizing 
reflection from the ground/air interface. Presently, 


Figure 1 - Schematic block diagram of a generic 
UWB GPR system 


The principal goal is to perform the needed 
analytical and design work to support the 
construction of an experimental/test facility, where 
targets can be buried in a sand medium and 
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illuminated by an ultra-wideband electromagnetic 
field, The scattered fields have to be observed and 
diagnosed for the purpose of building an electronic 
database containing the responses of a large number 
of buried objects. 


Generation of UWB EM fields 


For our UWB GPR system it was decided to use a 
solid state pulse generator with a double exponential 
voltage waveform. Such a waveform can be 
modeled by, 

V(t) = Volexp(- Br) - exp(- ar) }u(t), (1) 


with Va, a and B known constants, a >> B and 
u(t) the Heaviside step function. In our case the 
parameters are 
Vj= 9.635 kV, 2 
a=2. 10° sand B=25. 10° 5" (3) 


This gives rise to a Fourier spectrum, 


Bins 1 l 

Tj AP, a): el 
Since the magnitude of the far field is proportional 
to the derivative of the voltage waveform, a better 
radiated spectrum (enhanced bandwidth) is obtained 
when compared to sinusoidal mono-cycles. The 
pulser is designed to drive a 50 Q load, is fully 
protected agianst open and shortcircuit loads and has 
a pulse repetition frequency of 1 kHz. The pulse has 
a 10% to 90% rise time of about 100 ps and a full- 
width half-maximum (FWHM) time of 3 ns. 


Radiation of UWB EM fields 


Many commercially available wideband antennas 
have poor performance when it comes to radiating 
very short pulses. Part of the UWB GPR system 
design entailed the design of an IRA for ground 
penetrating radar applications. The upper frequency 
limit fupper of such an antenna is determined by the 
pulse rise time f,;;e while the lower frequency limit 
flower is determined by the dimensions of the 
reflecting dish. Pulse rise times of 25 ps and upper 
frequency limits of 20 GHz have been demonstrated 
in IRA designs. Experiments have shown that an 
antenna with a 4 m reflecting dish and a 60 kV 
switch with a 85 ps rise time has a bandwidth 
extending from 35 MHz to 4 GHz with a field 
stength of 4 kV/m at 300 m. The frequency range of 
interest in our case is from about 200 MHz to 3 


GHz. Hence the diameter of the reflecting dish was 
proposed to be in the order of 1 m. 


Soil chracterization 


Electromagnetic waves are damped considerably 
when traveling through a medium with a 
conductivity greater than zero. In our experimental 
set-up the angle of incidence and incident 
polarization are chosen such as to maximize the 
energy transfer from the transmitter into the ground. 
This is based on the technique of pulse matching 
into the ground developed in [10], for vertical polar- 
ization. If frequencies are much larger than the 
relaxation frequency of the ground (frerar = O/E), 
there is an angle of incidence at which there is very 
little reflection of the incident wave from the 
air/ground interface. This is called the high- 
frequency Brewster angle. For f >> frerax it can be 
shown that reflection is minimal (theoretically zero) 
when, 


tan(Ws) = Ve, (4) (5) 


When 5 = 0.01 S/m and e‚ = 10, then the relaxation 
frequency of the ground fret: = 113 MHz. Hence 
when the conductivity 5 < 0.01 S/m, which is quite 
typical, the radiated spectrum of the IRA will be 
transmitted into the ground with minimum 
reflection. The next table shows some examples of 
high-frequency Brewster angles at various values of. 
permittivity. 


e [2 4 6 8 10 
Wa | 54.73 6342° 67.19° 70.53° _ 72.45° 


Reception and synthesis of scattered fields 


The receiving system consist of: 

® __ receiving antenna, similar to transmitting 
antenna; 

® processing unit (data acquisition). 


The processing unit has to perform two tasks: 

* recording and preprocessing the incoming 
signal, reducing noise, clutter and distortion 
caused by the signal path; 

* perform some signal processing on the 
scattered field in order to be able to identify the 
buried objects. 


Data acquisition and processing is done by a 
sampling (digitizing) oscilloscope able to perform 
measurements on signals with bandwidth of up to 20 
GHz. The input signal cannot be measured reliably 
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at just one instance, which is why a sampling 
oscilloscope employs a repetitive sampling 
architecture. This in turn requires the transmitted 
pulse to be repetitive as well. Arithmetic averaging 
of the repetitive input signal can significantly reduce 
the noise (and jitter) on the received scattered field 
signal. Gating allows isolation of the part of the sig- 
nal due to the buried objects. 


Signal processing and identification 


Until 1980, almost all of the work in the field of 
groundpenetrating radar was directed towards the 
detection of buried objects. Only very few attempted 
the problem of target identification (or 
classification), which is a problem far more severe 
than the identification of aerospace targets by 
conventional radars where the target can literally be 
seen and the class of false targets is limited in scope. 
Underground there are varieties of undesired targets 
that complicate the task. Furthermore, the ground 
medium involved, is usually lossy, inhomogeneous 
and electrically weather-dependent. These problems, 
together with the presence of the air-ground 
interface, makes the task of subsurface target 
identification truly formidable. 

It was suggested by some authors that the so- 
called singularity expansion method (SEM) is 
one of the few methods capable of performing 
identification of objects [2]. This method is 
based on the theoretical observation that all 
objects have natural resonant frequencies that 
depend on their size, shape and the material 
decomposition only and not on the orientation 
of the object or on the direction of the incoming 
field. When a target is hit by a short 
electromagnetic pulse, the target will “ring” and ra- 
diate a scattered field that contains as prominent 
constituents a set of damped sinusoids. These 
sinusoids represent the natural resonant frequencies. 
The bandwidth of the ultra-wideband system must 
be large enough to cover all the resonant frequencies 
of interest. The SEM will be exploited in the current 
problem of classification of the buried object for the 
calculation of the internal resonances of the buried 
object by processing the reflections or 
backscattering from the objects. Some initial 
analysis with synthetic data has already been per- 
formed [1]. 


DESIGN AND CONSTRUCTION OF AN IRA 


The paraboloidal reflector antenna fed by a 
pyramidal horn has found wide-spread application in 


radar and communication engineering. The reflector 
antenna also has very useful characteristics when it 
is fed by two or four conductor transmission lines. A 
dispersion-less wideband antenna with a nearly flat 
radiating spectrum is desirable for short pulse 
applications. The reflector IRA employs a 
paraboloidal reflector fed by TEM lines [11, 12], 
and is an example of an aperture antenna. It is well 
known that the radiated field from an aperture 
antenna consists of a spatial integration of the 
aperture fields over the aperture, while the temporal 
behaviour of the aperture field is differentiated in 
the far field. In a practical situation, the illuminating 
field or the aperture field is a double exponential 
waveform and the radiated field then becomes 
impulse-like with a very large bandwidth ratio. 

The reflector IRA under consideration consists of 


odd, 
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eest Ce=rs 
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Figure 2 - An illustration of a reflector fed by a pair 
of coplanar conical TEM lines. 


a paraboloidal reflector fed by two pairs of coplanar 
feed plates as illustrated in figure 2. Coplanar feed 
plates are chosen over the more conventional facing- 
plate geometry to minimize the aperture blockage 
effects. To reduce the aperture blockage, the feed 
plates are required to be narrow, resulting in feed 
impedances of several hundreds of Ohms. The 400 
Q lines are connected in parallel resulting in a net 
feed impedance of 200 Q. The aperture area should 
be as large as practically possible, since the far field 
is proportional to the square root of this area for a 
constant voltage at the feed. The magnitude of the 
far field is proportional to the aperture area for a 
constant aperture field. The pulse generator has to be 
of the differential type to avoid common mode 
currents on the feed plates, which could distort the 
desired features in the far field. The driving voltage 
is V(t) = (Vo/2)u(t), where Vo is already.defined 
in eqn. 2. Since Eyar e« òV/òt , it; is desirable to 
maximize this rate of rise of the incident field or the 
voltage pulse. 
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Two identical reflectors were manufactured with the 
following characteristics. 


one piece, spun aluminum, paraboloidal surface 
diameter D =900 mm 

focal length F = 337.5 mm 

profile accuracy S 1.5 mm 

fu= FD = 0.375 


The desired rise time of 190 ps implies an upper 3 
dB frequency of 3.5 GHz which gives rise to 
wavelengths in air of 85.7 mm. The surface 
tolerance of the reflector is small compared to the 
shortest wavelength and hence it is acceptable. 


Next, we look at an estimation of boresight 
waveforms. For analysis purposes, one could 
consider a single (two conductor) coplanar feed, 
although in practice we used two such feed lines 
connected in parallel for a more uniform 
illumination of the reflector (see figure 2). When the 
reflector IRA was originally proposed [S], the 
boresight radiation was predicted to consist of a feed 
step followed by an impulse-like behaviour. It was 
also shown that the total area under these two parts 
of the radiated waveform (i.e. prepulse plus impulse) 
is zero. This means, that there is no DC component 
in the radiated waveform consisting of the prepulse 
and the impulse. This further implies that the portion 
of the radiated waveform after the impulse must 
have a net zero area in itself. The post impulse 
portion consists of diffracted signals from the feed 
plate and the circular rim of the paraboloidal 
reflector. A more recent analysis [6] has extended 
this result by chronologically considering the vari- 
ous temporal elements of the boresight radiation, 
which is illustrated in figure 3. Let us assume that 
the voltage pulse generator is switched on at t = 0, 
and the observer is at a distance r(= 2) to the right of 
the focal point of the paraboloid. These temporal 
elements are: 

© _Prepulse: feedstep E‚, (s 1) 

© Main pulse: impulse E‚» (r t) 

® _Postpulse: 

-_feedplate diffraction: E‚3(x t) (actually this 
consists of two parts originating form the 
plate edge (3a) and the plate itself (3b)) 

- edge diffraction from circular rim of the 
reflector, E‚4 (r‚ t) 

. Entire pulse constraints: low-frequency dipole 
moment radiation and no radiation at zero 
frequency (DC). 
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(ab) plate diffraction 


(3a) plate edge diffraction 


On-axis radiation from a canonical 
reflector IRA. 


Figure 3 


The far field E „(F,r) is given by, 
EDT EG), vee © 


E‚(rt)= Y E‚(r‚t) OD 


+low-frequency radiation from dipole 
moments resulting in time-integral 
constrains on entire pulse. 


Detailed, closed-form expressions for the 

components E,; (sf) in the above expression, 

are available in (7). The results of eqn. 6 is 

illustrated in figure 3. The time-integral 

constraints imposed by the low-frequency 

radiation are, 

i) the complete first-time integral of the radiated 
waveform must be zero, and 

ii) the second-time integral must be proportional 
to the late-time dipole moments. 


Feed plates, terminating and matching 


The feedplates form a parallel combination of 2 x 
400 Q conical TEM lines, giving the IRA an input 
impedance of 200 Q. Matching to the 50 Q pulse 
source is done by means of balun. Suitable high- 


„frequency, high-voltage, and high-impedance cables 


are difficult to find. Finally a 100 Q flexible cable 
with a PTFE dielectric was chosen. The latter is 
bound to introduce some dispersion, however the 
cable can be replaced by a better one in the future. 
Terminiation of each feed plate at the dish edge 
consists of a net DC impedance of 200 Q. The 
insulating spacer is a high density polyethylene slab 
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the dimensions of which are _optimised 
experimentally. Capacitance of the slab is trimmed 
by drilling holes (tuned with help of TDR 
measurements at the antenna terminal). The 
terminating resistor network, consist of a series- 
parallel network (5 by 5) of 1 W, 200 Q carbon 
composite resistors of 10% tolerance. Thus total 
heat dissipation is not a problem. Since the length of 
the resistor chain is > 10 cm, the voltage stand-off 
(in the order of 10 kV/cm) is not seen to be a 
problem either. 


Far field of the IRA 


If an aperture antenna of diameter D is illuminated 
by a CW field of frequency f, then the far field is 
determined by the distance r such that, 

2D° 
> 


is 


(8) 


which is obtained by requiring that A, < (À/16) 
where A, is the path difference of the edge ray and 
the central ray from aperture to observer. For pulsed 
antennas one can define a clear time ft, = A‚/c and 
require that f, < frise time Of the incident pulse. It is 
seen that 


Í zt DE 
hie d 3 (@) 


Requiring that fe < frise time » it is observed from 
table 1 that the observation point has to be 
greater than 8 m, to be in the far field, by a 
comparison of the clear time with the rise time 
of 100 ps. The radiated fields on the bore-sight 
are known in closed form in the far field [$, 6, 
7, 8, 9]. More recently, Russian researchers [14] 
have worked out near and far field expressions 
of the E-field on the bore-sight of an IRA. 
Using these expressions, we have estimated the 
on-axis E-field at various distances in the near 
and far field of the antenna. This is the subject 
of the next subsection. Although 8 m is the 
minimum distance to be in the far field, the 
actual experiments of illuminating targets can 
be performed at shorter distances in the near 
field. 

The pulse is broader and has more low-frequency 
components in the near field. In fact, very near the 
Table 1. Clear time f„ for a reflector diameter 
D=0.9 m 


r r+ A, A, 
(inm) (in m) 


te= Arle 


(in mm) (in ps) 


30 3.033 33 1100 
3:55 3529 29 96.6 
40 4,025 25 83.3 
50 5.020 20 66.6 
80 8.012 12 40.0 
10.0 10.010 10 33.3 


the antenna, for example at the feed point, the 
electric field is similar to the voltage pulse which is 
a double exponential waveform. The differentiation 
occurs in the far field and the transformation from 
near to far field is gradual and not abrupt. This 
transformation can be experimentally observed by 
measuring the electromagnetic field on-axis, as one 
moves away from the focal region. These values 


Table2. Near field and far field characteristics of 
the IRA (bore-sight) 


quantity near field far field 
distance R 5m 10 m 
cleartime 66 ps 33.3 ps 
prepulse -255 V/m -127 V/m 
prepulse duration 2.25 ns 2.25 ns 
impulse peak 5.9 kV 3.8 kV 
impulse duration 85ps 60 ps 


have been experimentally verified during initial 
measurements and compare favourably with the 
above table. Expressions for the radiated spectrum 
are available in literature as mentioned earlier. 
Calculations for our IRA results in the values 
tabulated in table 3. 


Table 3. Calculated bore-sight values for a 
distance R = 10 m 

E = 2.85. 10° V/m (+- 12%) 

fi GdB) = 80 MHz 

fu = 10 GHz 


BW-ratio = 125 

decades = 2.097 
Measurements indicate that due to practical 
restrictions the actual performance is not quite as 
good as theoretically expected. 
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UWB GPR EXPERIMENTAL FACILITY 


To be able to test the newly developed UWB 
GPR system and to perform experiments in a 
controlled environment, a full scale 
experimental/test facility was erected on the 
premises of TNO-FEL. This facility has been 
developed in cooperation with the Delft 
University of Technology and will be used to 
illuminate objects buried in a realistic sand 
medium with an ultra-wideband 
electromagnetic field, The electromagnetic field 
scattered by the buried object has to be recorded 
and diagnosed for the benefit of developing a 
computer database containing the responses of a 
large number of potential targets (ie. land- 
mines etc.). 


The experimental facility consists of a buried 
wooden box. The dimensions of this box are 10 m x 
10 m wide and 3 m deep. Special care has been 
taken not to use any metal parts in the construction 
of the box or in the vicinity thereof. The sand-box is 
filled with clean homogeneous river sand. In order 
to keep the condition of the sand in the box optimal 
and to prevent pollution from the outside (for 
example ground water) entering the box, a drainage 
system was installed and the inside of the box was 
covered with a watertight plastic lining. To prevent 
the weather from influencing the test conditions and 
to protect the measuring equipment, a large tent 
covers the entire site. While filling the box with 
sand, special care was taken to get a homogeneous 
profile. Current tests indicate that this was not 
entirely successful but should be good enough to 
work with. Later on we might choose to empty the 
box and refill it with a ground (soil) of a different 
composition. To facilitate the measurement of EM 
transmissions into the ground a square PVC tube 
running from the surface of one side to the bottom 
of the other side has been installed about one meter 
from the edge of the sand-box. 


Figure 4. The finished IRAs mounted on A-frames 
at the UWB GPR experimental facility. 
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IN MEMORIAM 


Op donderdag 28 januari 1999 overleed ons trouwe lid 
ir. P.D.C, Reefman op 66-jarige leeftijd als gevolg van 
een slopende ziekte. Kort na zijn afstuderen in 1963 te 
Delft als elektrotechnisch ingenieur trad hij in dienst bij 
het Dr. Neher Laboratorium van de PTT te 
Leidschendam, waaraan hij gedurende ongeveer tien jaar 
bleef verbonden. Tot de afsluiting van zijn actieve 
loopbaan ging hij daarna over naar de Hoofddirectie van 
PTT Telecom, met nadruk op zijn internationale werk. 


De opkomst van halfgeleidercomponenten, die als 
schakelmiddelen de functies van het magnetisch relais 
overnamen, had ingrijpende gevolgen voor wat toen werd 
aangeduid als de elektronische schakeltechniek. Ir. 
Reefman hield zich indringend bezig met de theoretische 
en praktische aspecten van die omschakeling, alsook met 
de gevolgen daarvan voor de beroepspraktijk van 
monteurs en technici. Hij maakte zich daarbij ook sterk 
voor standaardisatie van tekensymbolen, o.m. voor de 
verschillende poortschakelingen, trekkers, geheugens, 
vertragings-lijnen. De ontwikkeling daarvan vond in 
internationaal verband plaats en daarom had ir. Reefman 
zitting in een aantal internationale werkgroepen. 


Zoals uit het volgende blijkt toonde Ir. Reefman een 
grote betrokkenheid tot het onderwijs. 

Het verspreiden van de kennis omtrent de standaardisatie 
van schakelsymbolen, zowel in eigen 


kring als daar buiten, vroeg een bepaalde instelling die in 
hoge mate overeenkwam met het profiel van de docent. 
Daarnaast nam ir. Reefman gedurende vele jaren als 
examinator deel aan de schriftelijke en mondelinge 
examens van het NERG voor monteurs en technici. 
Ook schreef ir. Reefman op verzoek van de Vereniging 
Elektrotechnisch Vakonderwijs (VEV) leerstof voor op 
de PTT gerichte opleidingen in het kader van het 
leerlingstelsel. Zijn vaak herdrukt leerboek over de 
schakeltechniek werd op grote schaal ook daar buiten 
gebruikt. 

Sinds 1966 trad ir. Reefman op als gecommitteerde van 
het NERG voor het toezicht op het examen Hoger 
Elektronicus (later Elektronica Ingenieur) van Rens & 
Rens; vanaf 1974 tot vorig jaar als hoofdgecom- 
mitteerde. Van 1979 tot 1996 was hij ook als gecom- 
mitteerde van het NERG betrokken bij het examen 
Middelbaar Elektronicus van dit instituut. 


Wij kenden ir. Reefman als een beminnelijk mens. Hij 
wist steeds een goede persoonlijke sfeer te creëren en 
bracht zijn boodschap op een innemende en over- 
tuigende wijze. Zijn werk als gecommitteerde en 
examinator voor de Onderwijscommissie van het NERG 
genoot een hoge waardering. Zijn verscheiden betekent 
in vele opzichten een groot gemis voor het NERG. 


Wij wensen de familie Reefman sterkte bij de 
aanvaarding van dit grote verlies. 
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beurt moet wachten. Je kunt meteen je kennis en ideeën kwijt op grote projecten bij Signaal Hengelo of bij 
Signaal Delft. Omgekeerd kun je putten uit een enorme bron van kennis die wereldwijd gevoed wordt door 
meer dan 50.000 medewerkers van ons moederconcern Thomson-CSF. Met systemen als de Goalkeeper, de 
Flycatcher, APAR, Smart-L 5D radar en nog een aantal vindingen staat Signaal hoog in de octrooien-toptien. 

Maar op je lauweren rusten is er niet bij in de wereld van complexe radar- en gedistribueerde real-time 
gegevensverwerkende systemen, van command en control, sensoren en telecommunicatie. Daarom zoeken 
we (ervaren) elektrotechnici, technisch natuurkundigen en wiskundigen die als Radar Engineer en 
Systeemontwerper naam willen maken. LAAT HET MAAR ZIEN BIJ SIGNAAL 
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